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July  28  (Sat),  2012 

7:45  -  8:30  AM  Breakfast 

8:30  -  8:35  Welcome:  Hai  Wang  and  Yiguang  Ju 
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Chair:  Chung  K.  Law  (Princeton  University,  USA) 
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Zheng  Chen  (Peking  University,  China) 
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1 1 :05  —  11 :40 
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Session  3: 

13:30-14:05 

parameters: 

14:05-14:40 

14:40-15:15 

15:15-15:35 

15:35-17:30 


Predictive  Theoretical  Elementary  Reaction  Kinetics  and  its 

Role  in  Combustion  Modeling 

Stephen  J.  Klippenstein  (Argonne  National  Lab,  USA) 

Kinetic  Mechanism  Development  for  Hydrocarbons  and 

Oxygenated  Fuels 

Henry  Curran  (National  University  of  Ireland,  Ireland) 
Philippe  Dagaut  (CNRS  Orleans,  France) 

High  Pressure  Kinetic  Mechanisms  for  Hydrogen  and 

Hydrogen  Syngas 

Michael  P.  Burke  (Argonne  National  Lab,  USA) 

Frederick  L.  Dryer  (Princeton  University,  USA) 


Lunch 


Uncertainty  Quantification 

Chair:  Bill  Green  (IVIassachusetts  Institute  of  Technology,  USA) 
Determination  of  the  joint  uncertainty  domain  of  rate 

The  case  study  of  hydrogen  combustion 
TamasTuranyi  (Eotvos  University,  Hungary) 

The  role  of  uncertainty  Quantification  within  kinetic  model 

improvement  in  combustion 
Alison  S.  Tomlin  (University  of  Leeds,  UK) 

Mechanism  Optimization 

Hai  Wang  (University  of  Southern  California,  USA) 

Coffee  Break 

Group  Discussion  Session:  Identification  of  roadmap, 
challenges,  &  collaborations 
Elementary  Kinetics:  Chair,  Robert  S.  Tranter 


(Argonne  National  Lab,  USA) 
Mechanism:  Chair,  Matt  Oehlschlaeger 

(Rensselaer  Polytechnic  Institute,  USA) 
Flame  Chemistry:  Chairs,  Zheng  Chen 

(Peking  University,  China) 

Timothy  Ombrello 
(Air  Force  Research  Laboratory) 
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18:30-  Reception 


July  29  (Sun),  2012 

Session  4:  Combustion  Diagnostics 

Chair:  Nils  Hansen  (Sandia  National  Lab,  USA) 

8:30-9:05  Measurements  of  Elementary  Rate  Constants.  Ignition  Delay 
and  Species  Histories  in  Shock  Tubes 
Ronald  K.  Hanson  (Stanford  University,  USA) 
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research 
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Bielefeld,  Germany); 

Fei  Qi  (University  of  Chinese  Science  and  Technology,  China); 
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Elementary  Kinetics:  Chair,  Robert  S.  Tranter 

(Argonne  National  Lab,  USA) 
Mechanism:  Chairs,  Matt  Oehlschlaeger 
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(Peking  University,  China) 
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(Air  Force  Research  Laboratory) 

Presentations  of  Group  Discussion  by  Group  Chairs 


12:00-12:30 


Presentations  and  Workshop  Summary 


1.  Invited  Presentations  (11) 

2.  Workshop  summaries  (3) 

•  Elementary  Kinetics 

•  Flame  Chemistry 

•  Mechanisms 

3.  Workshop  presentation  (8) 


Laminar  Flames  and  the  Role  of 


Yiguang  Ju  and  Sang  Hee  Won 

Princeton  University,  USA 

Zheng  Chen 

Peking  University,  China 


1st  Flame  Chemistry  Workshop 


Advanced  Engines 

require  fuel  flexibility  and  work  near  kinetic  limit 


HCCI 


Synfuels  in  gas  turbine  en 


Low  temperature, 
High  pressure, 

Near  ign./ext .  limit, 
Multiple  fuels, 
Kinetic  limiting 


"Validated"mechanisms? 

•Hundreds  of  fuels 
•Different  structures 
elements  H,  C,  O,  N,  S...) 
•Thousands  species 
•Extreme  conditions 


Plasma  assisted  combustion 


Two  validation  targets 
•Homogeneous  ignition/reactor 


•Inhomogeneous  flames 


Flame  regimes  in  combustion 

How  does  chemistry  affect  flames? 


PIV 


Laser 


PLIF 


Turbulent  fuel  stream 
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Turbulent  flow  reactor 
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Reactive  flow 


•Thin  flame 
•Thickening  flame 
•Local  extinction 
•Re-ignition 


Propagating  edge  flame  in  a  mixing  layer 


Premixed  flame  front 


in  non-uniform  flow  field  with  complex  transport  and  chemistry  coupling 


>  "Flame 


Flames 

is  a  ignition/reaction  front  supported  by  thermal  and  species  transport 


Fuel  decomposition  ; 

•  H  abstraction  by  H 

•  Radical  termination] 


Fuel 

fragments 


Fuel 


Diffusion 

Fuel/Oxygen 
Radicals/Fragments, 
(e.g.  H  and  C2H4) 


Heat  release  rate 

OH  +  CO  =C02  +H 
HCO+OH=  C02+H20 


Non-uniform  species/temperature  distribution 


1.  Why  is  flame  chemistry  different  from  ignition? 


>Then,  what  is  the  role  of  transport  on  kinetics? 

•  Is  flame  chemistry  different  from  that  of 
homogeneous  ignition? 

•  How  does  transport  and  flame  chemistry  govern 
flame  extinction? 

•  How  does  transport  and  flame  chemistry  affect 
unsteady  flame  initiation  and  propagation? 

•  How  does  low  temperature  chemistry  change 
flame  regimes? 


Flames:  Different  fuels  have  different  extinction  limits 
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Tf=  500  K and  T0=  300  K 


aromatics 


O  A  0>  40  /-w 

O  A  <>0  A  CT" 

0  #  ft  ■ 

Ej'i4, 

*  ,«V  V  »  ** 


_  A  *♦ 
•#  ♦  w 


•• 

♦ 


How  to  decouple  chemistry  from  transport  and  fuel 

heating  value? 


0.05  0.1  0.15 

Fuel  mole  fraction  X, 


0.2 


Extinction  strain  rate  aE  [1/s] 


Ignition  vs.  flames:  n-alkanes  and  esters 
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Westbrook,  2010 
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Tf  =  500  K,  Tox  =  298  K 


O  Methyl  Formate 
O  Methyl  Ethanoate 
O  Methyl  Propanoate 

•  Methyl  Butanoate 

•  Methyl  Pentanoate 

•  Methyl  Hexanoate 
O  Methyl  Octanoate 
O  Methyl  Decanoate 


Kinetic  coupling  between  alkanes  and  aromatics 

Blending  toluene  into  n-decane:  Extinction  Limits 


Fuel  mole  fraction  Xf 


Won,  Sun,  Dooley,  Dryer,  and  Ju,  CF  2010 


Reaction  rate  [mole/cm3s] 


Kinetic  coupling  between  n-decane  and  toluene  in  diffusion  flames 
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60  %  n-decane  +  40  %  toluene 
nearextinction,Xf=  0.1,  a  =  1761/s 


C6H5CH2  +  H  =  C6H5CH3 


H  Diffusion  loss 


overall  decomposition 
of  n-dp£ane 


overall  decomposition 
of  toluene 


Cl  0H22=C2H5+C8H1 7-1 
Cl  0H22=nC3H7+C7H15-1 
Cl  0H22=CH3+C9H19-1 
C6H5CH3+H<=>C6H5CH2+H2 
C6H5CH3+0H<=>C6H5CH2+H20 
C6H5CH3+H<=>A1  +CM3 
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Won,  Sun,  Dooley,  Dryer,  and  Ju,  CF  2010 


Temperature  [K] 


Mass  burning  rate  (g/cmA2s) 


High  pressure  hydrogen  kinetics:  ignition  and  flames 
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-  Present  experiments 

■  Li  et  al.  (2007) 

■  Davis  et  al.  (2005) 

■  Sun  et  al.  (2007) 

■  Konnov (2008) 

O'Connaire  et  al.  (2004) 
Saxena  &  Williams  (2006) 
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Uncertainty  of  H02  related  kinetics  at  high  pressure. 


•Ignition  governed  more  by  chain  initiation  and  branching  rates 
•Flames  governed  more  by  branching  rate  and  heat  release  rate 
•Different  radical  pool  concentration  (H,  OH, ...) 

Bulke,  Marcos,  Dryer,  Ju,  CF,  2010 

Burke,  Chaos,  Ju,  Dryer,  Klippenstein,  IJCK  2011 


Difference  of  kinetics  in  homogeneous  reactor  and  flames 

Homogeneous  reactor  (800  K)  Fuel  consumption  by  radicals 
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OH  is  r  the  most  significant  radical  in  fuel  consumption 
Diffusion  flames  (~1600K)  puel  consumption  by  Radicals 


component 

%  by  uni- 
molecular 
decomposition 

%  consumed  by  radical  reaction 
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Methylbutanoate  [2] 
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1.  Won  etal.  CNF  159  (2012) 

2.  Dooley  et  al.  CNF  159  (2012)  1371-1384. 


2.  How  does  transport  and  flame  chemistry  govern  diffusion  flame 
extinction? 
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Diffusion  sensitivity 


Flames:  Different  fuels  have  different  burning  limits 
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How  to  decouple  chemistry  from  transport  and  fuel 

heating  value? 
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Fuel  mole  fraction  X, 
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A  generic  correlation  for  extinction  limit: 

Transport  weighted  Enthalpy  &  radical  index 

Theoretical  analysis  of  Extinction  Damkohler  number 
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Won  et  al.  CNF  159  (2012) 


A  General  Correlation  of  Hydrocarbon  Fuel  Extinction  vs. 

Transport  Weighted  Enthalpy  (TWE)  and  Radical  Index 
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•n-decane 
On-nonane 
On-heptane 
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•  n-propyl  benzene 
otoluene 

0 1, 2, 4-trimethly  benzene 
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R2  =  0.97 


Enabling  rapid  fuel 
screening! 


Tf=  500  K  and  To=300K 


0.5  1  1.5  2 

Rix[Fuel]xzlH^(Ml^ue/MI^/frose„)-1/2[cal/cm3] 


TWE  and  radical  index  for  predicting  extinction  limits  and 

synfuel  fuel  ranking  and  screening 


Radical  Index 
Representing  radical  pool 
High  temperature  reactivity 
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Ignition  Delay  vs.  Radical  Index  (real  fuel) 
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Extinction  of  diffusion  flame  in  counterflow  configuration 
Tf  =  500  K  and  Tair  =  300  K  @  1  atm 
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Dooley  etal.,  CNF  159  (2012) 


Consistent  in  high  temperature  reactivity 
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Measurements  from  RPI  (Oehlschlaeger),  Dooley  et  al.,  CNF  159  (2012),  Dooley  et  17 
al.,  CNF  (2012)  in  press 


Scaling  high  temperature  reactivity  of  methyl  esters 

Using  TWE  and  Radical  Index 


Dievart  et  al,  2012 

to  presented  on  Monday  at  34th  Symposium 


Tf  =  500  K,  Tox  =  298  K 

O  Methyl  Formate 
O  Methyl  Ethanoate 

O  Methyl  Propanoate  Different  heating  values 

•  v!c!!:y!  ®utfnoat®  Transport  properties 

•  Methyl  Pentanoate  ~  ~  ~ 

•  Methyl  Hexanoate 
O  Methyl  Octanoate 
O  Methyl  Decanoate 


Impact  of  alkyl  chain  length  on  methyl  ester  reactivity 

Methyl  Formate,  ROC  Methyl  Acetate,  R1C 


H  abstraction  reactions,  CH30C0  and 

CH30C(0)CH2  decomposition  reaction  to  presented  on  Monday  at  34th  Symposium 

rates:  large  discrepancies  (Xueliang,  2012) 


3.  How  does  transport  and  flame  chemistry 
affect  flame  initiation  and  propagation? 


Puzzle  of  high  altitude  relight: 
an  unresolved  ignition  problem  or  a  flame  problem? 


Altitude 

[i/p] 


Flow 

speed 


Flight  speed 


> 


Is  the  flame  speed  really  a  problem  for  relight? 


•  What  governs  the  ignition  &  Eig? 

•  What  are  the  chemistry  and 
transport  effects? 


_  Thermal  diffusivity 
Massdiffusivity 


oxygen 
Jet  fuel 


Eig  min:  Defined  by  flame  thickness,  5  (make  a  guess)? 


B.  Lewis  and  Von  Elbe  (1961),  Ronney,  2004,  Glassman  (2008) 

4  ,  11 

E  =-nd1,pC  (T  A-T  )oc  —  oc- 

ig  3  '  P V  ad  00  '  g  3 


3/2 


volume 


\ 

heat  capacity 


Larger  fuel  molecules  <=>  smaller  Ej 


*g 


•Ejg  mm*  Defined  by  stable  “flame  ball”  size? 

Zeldovich  et  al.  (1985),  Champion  et  al.  (1986) 

E,s=\^pCr(Tad-TJ~Le 

3  Larger  fuel  molecules  <=>  larger  Eig 


Theory:  Critical  Ignition  Size  vs.  Flame  Speed 


Assumptions  and  simplification: 

•  ID  quasi-steady  state.  Constant  properties 

•  One-step  chemistry 

•  Center  energy  deposition 


Tf- 

J  00 


r>-2  -U R  i  r>-2  -UL eR 

K  e  „  „_2  -u«  i  K  e 


Jr  2e  Vvdx 


-QR  e  =  — 

Le 


00 


=  exp 


|  r V““* 


Z  7>-l 


2  cr  +  (l-cr)r 


/ 


/? 


Ignition  by  heat  and  radical 
deposition  (qt=0.05) 


1:  qc  =  0.0 
2:  qc  =  0.5 
3:  q  =  0.675 
4:  qc  =  0.7 
5:  a„  =  0.73 


R 

n 


LeF  =  2.2 


Chen  et  al.  2011 
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Minimun  ignition  power,  Q  , 


Ignition  energy:  impacts  of  flame  chemistry  and  transport 


Cube  of  critical  flame  radius,  R3C  ^ 

O,  2 
</> 

Chen,  Burke,  Ju,  Proc.  Comb.  Inst.  Vol.33,  2010  i 

w  15 

15 

o 

'£  1 
'k. 

O 

0.5 


Fuel 

Mean 

Radical 

molecular 

Index 

weight 

JP8  POSF  6169 

153.9 

0.80 

SHELL  SPK 

POSF  5729 

136.7 

0.85 

@  1  atm 

Unburned  Temperature  =  450  K 
Fuel/Air  (21%  02)  mixture 


•  JP8  POSF  6169 

•  SHELL  SPK  POSF  5729 


0.6 


0.8 

Equivalence  ratio  ^ 


Won,  Santer,  Dryer,  Ju,  2012 


0.7 


0.9 
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Unsteady  flame  initiation 

Three  different  flame  regimes  (n-heptane/air) 


-  Regime  I 

•  Spark  assisted  ignition  kernel 

-  Regime  II 

•  Weak  flame  regime 

from  sparked  driven  ignition 

kernel  to  normal  flame 

-  Regime  III 

•  Self-sustained 
propagating  normal  flame 


Rapid  rise 
5.7  ms  2  ms 


Kim  et  al,  2012,  to  be  presented  on  Wednesday  at  34th  Symposium 


Flame  thickness  [mm] 


Rapid  change  of  flame  structure  in  flame  initiation  process 

(lean  n-heptane/air) 


0.006 

0.005 

0.004 

0.003 

0.002 

0.001 

0.000 

150  200  250  300  350  150  200  250  300  350 

Stretch  rate,  K  [sec'1]  Stretch  rate,  K  [sec'1] 

Completely  different  flame  structures! 

Kim  et  a  I,  2012,  to  be  presented  on  Wednesday  at  34th  Symposium 


Max  species  mole  fraction 


Normalized  flame  speed,  Sb  /  S 


Will  a  model  for  flame  speed  predict  the  unsteady  flame  initiation? 

An  example:  Lean  n-heptane  flame  initiation 


Critical  flame  initiation  radius  ( Rc )  >10  mm 


4.  How  does  low  temperature  flame  chemistry  affect 
flame  initiation  and  propagation,  and  stabilization? 


•Low  temperature  chemistry  (multi-stage  ignition) 
•Plasma  assisted  low  temperature  ignition,  extinction 


Location  of  maximum  heat  release  (cm) 


Multi  flame  regimes  in  HCCI  ignition  n-heptane:40  atm,  T=700  K 


1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

0.000  0.005  0.010  0.015 


Low  temperature  ignition 


Sf=25.6  m/s 


Sf=27.5  cm/s 
Low  temperature  flame  dominated 
double  flame  (decoupled) 

Sf=15.3  cm/s______ — *  Transition 

Single  high  temperature 
High  temperature  flame  flame  front 

dominated  double  flame  (coupled) 


LTI 

at  wall 


Hot  ignition 


Time  (s) 


Ju  et  al.  33rd  symposium  2011 


New  low  temperature  flame  regime  in 
Plasma  assisted  combustion 

(Counterflow  DME/02/He  ignition) 


emissions 

ch3o  +  oh  ->  ch2o*  +  h2o 

RCH20  +  OH  CH20*  +  ROH 

Sun  et  alv  Friday  34th  symposium  on  combustion,  2012  R:  organic  radicals 


OH  number  density  (cm" ) 


Kinetic  enhancement  of  plasma  assisted  ignition:  Change  of  S-curve 


He/02  =  0.66:0.34  and  0.38:0.62  ,  P  =  72  Torr,/=  24  kHz,  a  =  400  1/s 


7xl015 
6xl015 
5xl015 
4xl015 
3xl015 
2xl015 
lxlO15 

0.05  0.10  0.15  0.20  0.25  0.30  0.35 

Fuel  mole  fraction 


Plasma  assisted  LTC 
Plasma  generated 


Sun  et  al.,  34th  symposium  on  combustion,  2012 


Role  of  kinetics  on  PAC  at  low  temperature? 


Conclusion 


*  Flames  chemistry  differs  from  homogeneous  ignition  in  diffusion,  fuel 
decomposition,  radical  pool  production/consumption. 

•Low  temperature  and  unsteady  combustion  processes  lead  to  new 
different  flame  regimes  and  structures. 

•Flame  initiation  and  extinction  are  strongly  affected  by  both  transport 
and  chemical  kinetics. 

•Transport  weight  enthalpy  and  radical  index  are  developed  for 
predicting  extinction  limit  and  ranking  fuel  reactivity 

•Large  uncertainties  in  elementary  reaction  rates  of  kinetic 
mechanisms  for  simple  fuels  exist  in  extreme  conditions. 

•A  validated  mechanism  using  flame  speeds  fails  to  predict  unsteady 
flame  transition  and  the  critical  flame  radius. 
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Motivation 


Advanced  combustion  strategies  rely  on 

-  Low/moderate  temperature  combustion 

-  High-pressure  operating  conditions 

-  (Ultra)Lean  and  stratified  combustion 

-  Emerging  and  alternative  fuel  combustion 

Challenges 

-  Shift  from  mixing-controlled  to 
kinetics-controlled  combustion  regime 

•  Increasing  relevance  of  ignition- 
kinetics  and  low-temperature 
chain-branching  reactions 

•  Increasing  significance  of  turbulence 
and  turbulence/chemistry  interaction 

•  Finite-rate  chemistry  effects 

-  Operation  near  stability  limit 


---125mm 

- 35mm 
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Motivation 


•  Objective 

-  Development  of  high-fidelity  combustion  for  prediction  of 
turbulent  reacting  flows  under  consideration  of 

•  Finite-rate  chemistry 

•  Turbulence/chemistry  coupling 

•  Transient  combustion-dynamical  processes 

•  Relevance 

-  Identify  and  isolate  combustion-physical  processes 

-  Combustor-design,  control,  and  optimization 

-  Guide  experimental  instrumentation 
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Overview 


•  Motivation 

•  LES-combustion  modeling 

-  Flamelet-based  formulation 

•  Part  1:  Modeling  and  simulation  of  combustion-physical 
processes:  LES  of  lifted  vitiated  flames 

•  Part  2:  Guide  experimental  instrumentation:  Turbulent 
inhomogeneities  and  facility-effects? 

•  Summary  and  conclusions 
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Overview 


•  LES-combustion  modeling 
-  Flamelet-based  formulation 
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LES  Combustion  Modeling 


•  LES  Flamelet-based  combustion  models 


-  Representation  of  turbulent  flame  as 
unsteady  reaction-diffusion  layer  that  is 
embedded  in  turbulent  flame 

-  Interaction  of  flame  structure  with 
turbulent  environment  leads  to 
stretching,  deformation,  and  extinction 
of  flame 
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LES  Combustion  Modeling 


•  LES  flamelet-based  combustion  model 


2500 
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LES  Combustion  Modeling 


•  Flamelet-structure  in  turbulent  reacting  flows 

-  Analysis  Tools:  DNS-database1  of  reacting  jet-in-cross-flow 

•  Fuel:  N2-diluted  H2-jet,  350  K 

•  Oxidizer:  Air,  750  K 


•  Extract  instantaneous  local  flamelet  structure  from 
DNS-database 


1  Grout  et  al.  Proc.  Combust.  Inst.  2011;  Kolia  et  al.  US  Nat.  Meeting  2011 
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LES  Combustion  Modeling 


•  Evolution  of  lD-flamelet-elements  in  JICF 


LES  Combustion  Modeling 


•  LES  Flamelet-based  combustion  models 


-  Parameterization  of  combustion  process 
in  terms  of  reduces  set  of  scalars 

-  Account  for  detailed  chemistry 

-  Tabulation  of  reaction  chemistry 

-  Consideration  of  turbulence 
chemistry  coupling 
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Overview 


•  Part  1:  Modeling  and  simulation  of  combustion-physical 
processes:  LES  of  lifted  vitiated  flames 
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Autoignition  in  Turbulent  Flames 


Modeling  challenges  in  predicting  autoignition  in 
turbulent  flames 

-  Autoignition  is  transient  process; 
requires  accurate  description  of 
temporal  flame-evolution 

-  Flame  stability  and  ignition 
dynamics  strongly  dependent  on 
scalar  mixing  and  flame/turbulence 
interaction 


Modeling  approach1 


-  Autoignition  requires  consideration  of  transient  species 
formation, 

described  by  unsteady  flamelet  equations 

-  Turbulence/chemistry  interaction:  Presumed  PDF-closure  to 
consider  effects  of  subgrid-mixing  and  unresolved  flame 
structure 


1  Ihme  &  See,  Combustion  &  Flame,  157,  2010 
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Autoignition  in  Turbulent  Flames 


•  Conditions  for  flame-ignition  in  diffusion  flames 

-  Autoignition  is  transient  process 

-  Sufficiently  low  scalar  dissipation  rate 

-  Flame  ignition  occurs  under  conditions  corresponding  to 
“most-reactive  mixture” 

-  Build-up  of  radical  pool  through  chain-branching  reaction 
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Autoignition  in  Turbulent  Flames 


•  Experimental  configuration 

-  Lifted  flame  in  vitiated  co-flow 

-  Fuel:  methane/air  1:2 

-  Co-flow  temperature:  1350  K 

-  Co-flow  composition  from  premixed 
H2-Air  reaction  product 

•  Computational  setup 

-  Grid:  2.5  Mio  grid  points 

-  Reaction  Chem.:  GRI  2.11, 

(also  used  GRI  3.0,  USC-mech  II) 

-  5-dimensional  chemistry  table  with 
grid-refinement 


1  Cabra,  Chen,  Dibble,  Karpetis,  Barlow,  CF,  143,  2005 
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Autoignition  in  Turbulent  Flames 


8 


6 

^4 

CL 

2 

0 


•  Ignition  conditions:  low-strain  region  at  most- 
reactive  mixture  composition 

•  Ignition  occurs  primarily  in  diffusion  regime 

•  Location  of  flame-base  controlled  by  H02- 
radical  pool  that  is  formed  upstream  of  flame 


Ignition  Mode  at  Flame  Base 


-1  -0.5  0  0.5  1 

Flame  Index,  FI 
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Autoignition  in  Turbulent  Flames 


•  Effects  of  turbulence  and  scalar  mixing 
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Autoignition  in  Turbulent  Flames 


Instantaneous  temperature  field 
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Autoignition  in  Turbulent  Flames 
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Autoignition  in  Turbulent  Flames 


•  Centerline  profiles 
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Autoignition  in  Turbulent  Flames 


Radial  profiles 
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Summary  and  Conclusions 


•  LES-modeling  of  lifted  vitiated  flames 

•  Key  modeling  components 

-  Transient  flame  evolution 

-  Accurate  description  of  turbulent  mixing  and  scalar  dissipation 
rate 

•  Combustion-physical  insights 

-  Transient  flame  evolution 

-  Identified  significance  of  flame/turbulence  interaction 

->  Homogeneous  reactor-model  under-predicts  ignition  onset 
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Overview 


•  Part  2:  Guide  experimental  instrumentation:  Presence  of 
turbulence  in  chemical-kinetics  facilities? 
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Turbulent  Inhomogeneities  and  Facility  Effects 


•  Question:  Can  we  apply  "lessons-learned"  from  LES 
simulations  to  characterize  experimental  facilities? 

-  Shock-tubes 

-  Flow-reactors 

-  Rapid  compression  machines 


•  Source  of  non-idealities  in  experimental  facilities1'2-3 


1  Petersen  et  al.,  Comb.  &  Flame,  149,  244  (2007);  2  Chaos  &  Dryer,  Combust.  Sci.  and 

Tech.,  180,  1053  (2008);  3  Burke,  Chaos,  Dryer,  &  Ju,  Comb.  &  Flame,  157,  618  (2010) 
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Turbulent  Inhomogeneities  and  Facility  Effects 


•  Research  Objectives 

-  Use  high-fidelity  simulation  and  non-equilibrium 
formulation  to  isolate  parametric  contributions  of  non¬ 
idealities  in  experimental  facilities 

-  Research  emphasis 

•  Identify  parametric  sensitivities 

•  Reconcile  observed  differences  between  experiments  and 
detailed  model-formulations 

-  Facilities 

•  Shock-tube 

•  Flow-reactor 

•  Rapid  compression  machine 
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Turbulent  Inhomogeneities  and  Facility  Effects 

Shock-Tube 


•  Non-ideal  processes  in  shock-tubes 

1.  Non-ideal  rapture  of  diaphragm* 

•  Finite  opening  time  of  diaphragm 

•  Contribution  to  shock  attenuation:  30% 

2.  Boundary  layer  growth $ 

•  Formation  of  viscous  boundary  layer 
behind  initial  shock 


3.  Shock  reflection  and  bifurcation0/o 

•  Lift-off  of  boundary  layer  resulting 
in  formation  of  separation  region 


Strong  ignition/ 


4.  Inhomogeneous  ignition  and 

weak-to-strong  ignition  transition $ 

•  Ignition  proceeds  as  multi-dimensional 
heterogeneous  process 


#  Petersen,  E.  L.  &  Hanson,  R.,  Shock  Waves,  10  (2001);  $  Mirels,  H.,  NACA-TN  3278  (1956); 
Mirels,  H.,  NACA-TN  3401  (1955);  %Yoo,  Mitchell,  Davidson,  Hanson,  Shock  Waves,  21,  2011 
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Turbulent  Inhomogeneities  and  Facility  Effects 

Shock-Tube 


•  Modeling  challenges  in  simulating  shock-tubes 

-  Disparity  of  spatial  and  temporal  scales 


150  Re2J7  L  ~ 

i _ —L _ i 

i 

Geometric  shock- 
tube  domain 


ReT  M 
Mt 
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l 


I 

Shock 

Thickness 
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•  Solution  method:  Adaptive  mesh 
refinement 

-  AMR  exploits  multiscale  nature  of 
hydro-dynamic  problem  by  locally 
adjusting  computational  effort  to 
maintain  uniform  level  of  accuracy#'* 


#  Deiterding,  R.,  Computers  &  Structures,  87,  769,  2009;  $  Berger,  M.  &  Colella,  P.,  JPC,  1982 
(1988);  %  C.  Pantano,  et  al.  JCP,  221  (1)  (2006)  63-87. 
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Turbulent  Inhomogeneities  and  Facility  Effects 

Shock-Tube 

•  Shock-bifurcation 

-  Simulation  of  Ar-diluted  H2/02  mixture  at  5  and  10  bar 
pressure 

^Relevant  condition  for  weak  and  strong  ignition  regime 
-> Adiabatic  and  isothermal  wall  conditions 

-  Shock  tube  setup 

•  Length:  lm 

•  Diameter:  5  cm 

•  Helium  in  driver  section 

-  Target  condition: 

•  T5=1100  K,  p5=10  bar 

-  Chemical  mechanism: 

•  Burke  et  al.1  (2011) 


25.0 

24.0 

23.0 

1  22.0- 1 
o 

J  21 . 0 
•3  20.0-1 

19.0  A 
18.0 
17.0  A 


Simulation 


Experiment2 


1  Burke,  Chaos,  Ju,  Dryer,  Klippenstein,  IJCK,  44,  2011 

2  Yoo,  Mitchell,  Davidson,  Hanson,  Shock  Waves,  21,  2011 
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Results:  Shock  Bifurcation 


•  Shock-bifurcation 


-  Instantaneous  temperature  evolution  shows  rich  flow-field 
structure:  Boundary  layer  separation,  Shock-bifurcation, 


JT-Axis 


Toluene  PLIF 
Measurements# 


#  Yoo,  Mitchell,  Davidson,  Hanson,  Shock  Waves,  21,  2011 
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Turbulent  Inhomogeneities  and  Facility  Effects 

Shock-Tube 


•  Ignition 

-  Isothermal  wall 


Var:  Temperature 
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-  Observations: 

•  Ignition  is  initiated  at  end-wall 

•  Flame  propagation  towards  unburned  mixture  (region  of  favorable 
pressure  gradient) 

->  non-homogeneous  ignition 
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Turbulent  Inhomogeneities  and  Facility  Effects 

Flow-Reactors 
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Peschke  &  Spadaccini,  (1985);  Santoro  (2009) 
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Turbulent  Inhomogeneities  and  Facility  Effects 

Flow-Reactors 
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Peschke  &  Spadaccini,  (1985);  Santoro  (2009) 
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Turbulent  Inhomogeneities  and  Facility  Effects 

Flow-Reactors 


Experimentally  observed  stochastic 
ignition  suggests  sensitivity  to 
initial  conditions 

-  Mixture  composition 

-  Temperature 

•  Unsteady  heating 

•  Wall-heat  losses 

•  Temp-difference  btw.  fuel  and  oxidizer 

Consider  inhomogeneities 


-  Equivalence  ratio:  sample  from 
experimentally  determined 
beta-distribution 

-  Temperature  fluctuations:  Sample 
from  Gaussian  with  specified  V 

Use  fully-developed  turbulent 
pipe-flow  at  Re  =  104 


1  Santoro  (2009) 

2  Samuelsen  et  al.  (2003) 
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Turbulent  Inhomogeneities  and  Facility  Effects 

Flow-Reactors 


•  Mixture  variation:  4>=0.4;  4>  -0.135 

•  Temperature  variation:  T=850  K;  T'  =  {0,  25,  50,  75}  K 


T’  =  0K 


T’  =  25  K 


0  0.25  0.5  0.75  1  1.25 
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Turbulent  Inhomogeneities  and  Facility  Effects 

Flow-Reactors 


•  Mixture  variation:  4>=0.4;  4>  -0.135 

•  Temperature  variation:  T=850  K;  T'  =  {0,  25,  50,  75}  K 
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Summary  and  Conclusions 


•  Turbulence/chemistry  coupling  processes 

-  Increased  relevance  for  low-Damkoehler/high-Karlovitz  combustion 
processes:  oxygen-diluted  comb.;  autoignition;  preheat-comb. 

-  Turbulence  promotes  mixing,  exchange  of  radicals  and  enthalpy 

-  Ignition  occurs  at  preferred  sites:  "most-reactive"  mixture  and  regions  of 
low  strain 

•  Validated  high-fidelity  LES  combustion  models  have  been 
developed  and  are  available  to  accurately  capture  ignition 
processes 

-  Models  rely  on  experimental  data 

•  Simulations  can  assist  and  complement  experimental 
investigations 

-  Identify  experimental  sensitivities 

-  Guide  potential  modifications  to  mitigate  facility  effects 

-  Reconcile  discrepancies  btw.  experiments  and  theory 

-  Example:  Turbulence/chemistry  coupling  in  shock-tubes  and  flow-reactors 
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Micro  flow  reactor  with  prescribed  temperature  profile 

Toward  fuel  Indexing  and  kinetics  study 
based  on  multiple  weak  flames 
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K.  Saruwatari,  S.  Suzuki,  T.  Kamada,  X.  Li,  Y.  Kizaki,  T.  Onishi  and  H.  Takahashi 
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Micro  flow  reactor  with  prescribed  temperature  profile 
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Micro  flow  reactor  with  prescribed 

temperature  profile 
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Flame  behaviors  in  a  micro  flow  reactor 
with  a  prescribed  temperature  profile 
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Three  kinds  of  flame  responses 
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Triple  weak  flames,  n-heptane 
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LTO:  CH20,  H202,  CO,  CH4  produced 
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Weak  flames  at  different  RON 
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Hori  etal.,  CNF159  (2012) 
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Weak  flames  at  various  RON 


Data  collected  at  slightly  different  conditions, 
i.e.,  flow  velocity,  exposure  time,  temperature  profile 
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Conclusions  and  future,  l  of  2 


Micro  flow  reactor  with  prescribed  temperature  profile  was  introduced 

-Three  kinds  of  flame  response  (S-shaped) 

-Weak  flame  corresponds  to  ignition  branch 
-Multiple  weak  flame  utilized  for  fuel  characterization 
-Multiple  weak  flames  at  elevated  pressures 
-Appearances  of  weak  flame  correlated  with  RON 
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Conclusions  and  future,  2  of  2 

Diesel  fuels  and  Cetane  numbers 

n-cetane  (hexadecane,  C16H34) 

/'so-cetane  (2,2,4,4,6,8,8-heptamethylnonane,  C16H34) 
/i-decane  (C10H22) 

a-methylnaphthalene  (C^H^)  — ►  (Suzuki  et  al.,  5E06) 

Ethanol  —*■  (Nakamura  etal.,  1E02) 

Syngas 

Oxyfuel  combustion  — ►  (Li  et  al.,  WIPP) 

Effect  of  surface  reaction  — ►  (Kizaki  et  al,  WIPP) 

PAH  and  soot 

Optical  diagnosis 

LIF,  CRDS  for  precise  species  profile  measurements 

Higher  pressures 

High  pressure  chamber  (  up  to  20  bar)  fabricated 
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Micro  flow  reactor  now  commercially  available 
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1.  Theoretical  Studies  Motivated  by  Modeling 

•  CO  +  H02,  H  +  H02,  CH3OH  +  H02 

•  NNH  Chemistry 

•  CH2  +  02 

•  c3h6ooh  +  o2 

•  C4H9OH,  C4H9O,  C4H9O  +  02 


Uncertainties 

•  Electronic  Energies 

•  Vibrational  Frequencies 

•  Torsional  Modes 

•  Energy  Transfer 

•  Coupling  of  TST  and  Energy  Transfer 
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Validated  Against  Ignition  Delays,  Flame  Speeds,  Species  Profiles 


ch3oh  +  ho2  =  ch2oh  +  h2o2 


H202(+M)  =  20H(+M) 


CH3OH  +  OH  =  CH3O  +  H20  • 


SH3OH  +  OH  =  CH2OH  +  H20 


|  H  +  02  =  O  +  OH 
|  HCO  +  M  =  H  +  CO  +  M  - 


OH  +  CH3(+M)  =  CH3OH(+M) 


|  CO  +  OH  =  C02  +  H  | 
CH2OH  +  M  =  ch2o  +  H  +  M  - 


ch3oh  +  h  =  ch2oh  +  h2 


HCO  -f  P2  =  CO  +  HQ2 
CH20  +  OH  =  HCO  +  H20 


ho2  +  oh  =  h2o  »  o2  - 


2H02  =  H202  +  02 


-0.8 


M  VPFFT1000 
□  Ignition  delay 
Flame  speed 


-0.4  0  0.4 

Sensitivity  coefficient 


0.8 


How  to  Improve  for  Engine  Conditions  -  High 
Pressures/Non-Dilute? 


Global  uncertainty 
analysis 

•  Indicates  reaction  with 
dominant  contribution 
to  overall  uncertainty 

•  Theoretical  predictions 

of  selected  rate 
coefficients  & 

•  How  accurate? 

Klippenstein,  Harding, 
Davis,  Tomlin,  Skodje 
Proc.  Comb.  Inst.  33, 

351  (2011). 


Variance  Analysis  for 
the  Ignition  Delay 

0.8 
0.7 
0.6 
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0.3 
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0.1 
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0  20  40  60  80 
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CH  OH+O  =CH_OH+HO_ 
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HO.+HO  =Ho0.+0_ 
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H202(+M)=OH+OH(+M) 


h. 
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100 


Transition  State  Theory 

Complete  torsional 
sampling  to  find 
torsional  minima 

Torsional  modes  treated 
as  uncoupled  1- 
dimensional  hindered 
rotors 

Asymmetric  Eckart 
Tunneling 

Other  modes  - 
harmonic  oscillators 
and  rigid  rotors 


CH3OH  +  ho2 


CO 


Change  in  ignition  characteristics 

Predicted  Ignition  Delays 


Reaction  Index  i 


Initial  0.11 
Final  0.40 


Initial  2.0 
Final  0.70 


How  Accurate  are  our  Predictions? 


■  Currently  a  Factor  of  2-3  Uncertainty 

■  Compare  with  Experiments 

-  Experience  from  many  comparisons 

■  Analysis  of  uncertainty  for  underlying  factors  in 
theoretical  analysis 

-  Barrier  Height  Uncertainties 

-  Vibrational  Frequency  Uncertainties 

-  Anharmonicity  Corrections 

■  Do  Mechanisms  Improve? 

■  Trying  to  reduce  uncertainties  to  about  1.2-1. 3 


Reactions  with  Barriers 


Boltzmann  Factors  for 
Given  Errors  (kcal/mol) 


Combustion  Thermochemistry  Database 

HEAT  (Stanton  and  coworkers);  Wn  (Martin  and 
coworkers);  Focal  Point  (Allen  and  coworkers) 

CCSD(T)/TZ  optimizations  and  Frequencies 

CCSD(T)/CBS  from  CCSD(T)/AQZ’,A5Z’ 

CCSDT(Q)/DZ 

Core-Valence  CCSD(T)/CBS;  TZ  and  QZ 
Relativistic  from  DKH  with  CI/TZ 
DBOC  from  HF/cc-pVTZ 

Anharmonic  corrections  from  B3LYP/6-31 1 ++G** 

Heats  of  formation  relative  to  H2,  CH4,  H20,  NH3 
5  Heavy  Atoms 


1 

Size  of  Individual  Corrections 


Mean 

MAD 

RMSD 

Variance 

CCSDT(Q) 

-0.36 

0.42 

0.62 

0.46 

Core-Valence 

0.54 

0.54 

0.60 

0.28 

Anharmonicity 

0.47 

0.49 

0.59 

0.32 

Relativistic 

-0.16 

0.16 

0.20 

0.11 

DBOC 

0.10 

0.11 

0.14 

0.09 

Number 


Accuracy  vs  ATcT  from  Ruscic 
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Heat  of  Formation  Errors 


I  1  1  1  I  1  1  1  I  1  1  1  I  1  1  1  I  1  1  1  I 


0 


With  Corrections 
Without  Corrections 


0.05  0.15  0.25  0.35  0.45  0.55  0.65  0.75 


kcal/mol 


H,  O,  OH,  02,  H02,  H202, 
03,  C,  CH,  CH2,  1CH2,  CH3, 
CO,  HCO,  COH,  CH20, 
HCOH,  CH20H,  CH30,  C02, 
HOCO,  HC02,  OCHOH, 
CH302,  CH300H,  C2, 

C2H2,  H2CC,  C2H3,  CCH3, 
C2H4,  CHCH3,  C2H5,  C2H6, 
CCO,  1CCO,  HCCO, 

CH2CO,  HCCOH,  CH3CO, 
CH3CHO,  CH2CHOH, 
CH3CHOH,  C2H40H, 
CH3CH20,  C2H50H, 
CH30CH3,  OCHCHO, 
CH3C(0)0H,  CH2CCH, 
CH3CCH,  CH2CCH2,  - 
CH2CHCH-,  CH3CHCH2, 
CH3CH2CH3,  N,  NNH, 
NHNH,  H2NN,  NHNH2, 
NH2NH2,  NO,  HNO,  NOH, 
H2NO,  HNOH,  H2NOH,  N20, 
N02,  N03,  CN,  HCN,  CNH, 
H2CN,  HCNH,  CH3N, 
CH2NH,  CH3NH,  CH2NH2, 
CH3NH2,  NCO,  HNCO, 
HCNO,  NCOH 


QCI//CC/TZ 


Rovibrational  Properties  Reactions  on  C3H7  PES 


QCISD(T)/CBS//B3LYP/6-31 1  ++G(d,p) 


- n->i 

- n->H+C3H6 

— 1 — 1 — 

- i->n 

— H+C3H6->n 

1 

. n->CH3+C2H4 

—  *■  -  i->C3H6+H 

• 

- CH3+C2H4->n 

-  -  a-  -  H+C3H6->i 

400  600  800  1 000 1 200 1 400 1 600 1 800  2000 


QCISD(T)/CBS//M06-2X//6-31 1  ++G(d,p) 


Temperature  (K) 


Temperature  (K) 


QCI//CC/TZ 


Rovibrational  Properties 


QCISD(T)/CBS//PT2/TZ 


QCISD(T)/CBS//CCSD(T)/ADZ 
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Temperature  (K) 


Temperature  (K) 


Coupling  of  Uncertainties 


AE  =  ±1 .0 
kcal/mol 

Av  =  1.5  factor 

v*  =  (50,850) 


CH3CH(OH)CH3  +  OH  —  CH3C(OH)CH3  +  h2o 


cm 


-1 


J.  Prager,  H.  N. 
Najm,  J.  Zador 
Proc.  Comb.  Inst. 
In  press  (2012) 
Paper  5D03 


SJK 


Direct  Sampling  of  Torsions 


CH  OH  +  HO  ->  CH  OH  +  H  O 


A 


k(s  ) 


Pressure  Dependence:  Comparison  with  Expt 


Propane  Decomposition  Acetaldehyde  Decomposition 


1000  K/T 

<AEdown>  =  100  (T/300)0-85  cm-1 


1000  /T(K) 

<AEdown>  =  1 50  (T/300)0-85  cm-1 


CoH 


2'  '3 


->  C2H2  +  H 


Energy  Transfer 
Jasper  &  Oana 


A  Priori  Predictions  of  Pressure  Dependence 


Trajectory  Simulations  of  Energy  Transfer 
A  Priori  P(E’->  E) 

ID  Master  Equation  Assumes  Strong  Collider  in  J  -  Not  Valid 
P(E’,J’->E,J)  AEdown  AJdown 
2-Dimensional  Master  Equation  (E,J) 


Coupling  of  Uncertainties:  n-propyl  +  0 


02+I1-C3H7 

ro2 
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N 
♦  «■» 

0.06 

OH+propene-oxide 

[R7] 

O 

0.04 

QOOH-1 

H02+propene 

[R8] 

c 

0.02 

7 — 

OH+propene-oxide 

[R9] 

0.00 

0.1  1  10 
normalized  rate  constant 

C.  F.  Goldsmith,  A.  S.  Tomlin,  S.  J.  Klippenstein,  Proc. 

Comb.  Inst.  34,  in  press  (2012)  Paper  2D09 
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Future  HC  fuels  -  many  sources 


C; 


>  Some  petroleum  will  still  be  available 

>  Oil  sands,  oil  shale 

>  Coal-to-liquids 

>  Fischer -Tropsch 

>  Natural  gas 

>  Hydrogen 

>  Bio-derived  fuels 

>  Ethanol,  butanol,  algae 


a  ^  ^  ^  ^ b  ^  cF 

O  O  O  Cp  O  GOO  Cj  0  GOO  G  G  G  G  G 


HYDRO WR EON  HUMBER 


>  Biodiesel  from  vegetable  and  animal  oils 


>  Chemical  kinetics  to  understand  and  simulate 


>  complex  behaviour  (ignition,  NTC,  cool  flames...) 
>  reactivity  (extent  of  conversion,  heat  release) 

>  product  /  pollutant  formation 
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How  well  an  oxygenated  fuel  works  depends  on  its 

molecular  structure 
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Reference  fuel  mixed  with: 
-■-MB 
-•-  DMM 
-a-  TPGME 
▼  Ethanol 
Methanol 
DME 

Dibuthylmaleate 


0.5 


LLNL 


Fraction  of  oxygen  in  fuel  by  mass 

Miyamoto  etal.  Paper  No.  SAE  980506  (1998). 

Westbrook  etal.  J.  Phys.  Chem.  A  (2006)  770:6912-6922.  ff j)| 
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Oxygenated  fuels 


>  Alcohols  (methanol,  ethanol,  propanol,  butanol) 

>  Ethers  (DME,  DEE,  EME,  MTBE,  ETBE) 

>  Esters  (methyl  and  ethyl  esters) 

>  Ketones  (acetone,  EMK,  DEK) 

>  Furans  (methyl  furan,  di-methyl  furan) 
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General  reaction  scheme 


C; 


High  temperature _ 


( 


6® 


.0° 


icctf 


p* 


p*  x  uel 


i 


„  ,  ,  p*scrs5Ton  .  +  02  _  • 

smaller  molecules  R  i  >  RO- 


olefins, 
cyclic  ethers,  <  i 
(3-scission 

Intermediate  temperature 


\7 


QOOH 


Low  temperature 

+  o2  chain 
branching 
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Fuel  decomposition  reactions 
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Ignition  delay  time  (ns) 


Propanol  isomers:  Comparison  of  reactivity 


0.5%  fuel,  2.25%  02  <|>  =  1 .0,  P  =  1  atm 

1000 


100 


5.0  5.5  6.0  6.5  7.0  7.5 

104K  IT 


■  iso  -propanol 
•  n  -propanol 


*  • 
um  «• 

ch3  ; 

CH3— CH — OH  ^  • 

*  CH3 - CH2— CH2— OH 


Johnson  ef  a/.  Energy  &  Fuels  (2009)  23  5886-5898. 
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Alcohol  molecular  elimination 


H— C— H 
H  -  C-OH 


J 

H 


H 


n-propanol  propene  +  water 

k  =  3.52  x  1013  0.  67300.  S1 


H, 

H—C-H 


H3C  -  C-OH 

i 

H 


H 


/'so-propanol  propene  +  water 

k  =  2.11  x  10140.  67300.  S'1 


Tsang,  W.  Int.  J.  Chem.  Kinet.  (1976)  8: 173-192. 
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Water  elimination  is  much  more  important  for  iso¬ 
propanol 


C; 


0.5%  /'so-propanol 
c|>  =  1.0,  T=  1600  K 
30%  fuel  consumed 


(5.6%) 


CH3 

H  (10%)  I  H  (7.7%) 

- -  H^C— CH— OH  - ► 

OH  (5.2%)  3  OH  (5  8%) 

0(1.6%)  OC3H7OH)  OH(58/o) 


163%) 


h3c—  ch=ch2  +  h2o 
(C3H6) 


Johnson  et  at.  Energy  &  Fuels  (2009)  23: 5886-5898. 
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General  reaction  scheme 


High  temperature 

Fuel 


olefins, 


cyclic  ethers,  <  i  QOOH 
(3-scission 


Low  temperature 


chain 

branching 


Intermediate  temperature 
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Sub-mechanism 
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||  OE  Gaillimh 


Reactivity  of  ethers 

•  1%  H2, 1%  02  in  Ar,  p5  =  1 .4-2.6  atm 


Yasunaga  etal.  Comb.  Flame  (2011)  758/1032-1036. 
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Ether  molecular  elimination 


C3 


H, 

H-C-H 

h3c-  C-O-CH3 

l 

H3C 
MTBE 

k  =  1.70  x  1014  0.  60800.  S'1 


+  H3C-OH 


isobutene  +  methanol 


H, 

H-C-H 

I 

H  qm-'C  O  OH2  CH3 

h3c 

ETBE 


+  H3C-CH2— OH 


isobutene  +  ethanol 


k=  1.70  x  10140.  60600.  s1 


Yasunaga  etal.  Comb.  Flame  (2011)  755/1032-1036 

MP4/cc-pVTZ//MP2/cc-pVTZ  level  of  theory  with  zero  point  corrections 
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Ether  molecular  elimination 


H 

1  ,H 

h3c  h-c 

h3c- c-°-ch2 

h3c 

ETBE 

y 

h^c-h 

I 

H2C-0-CH3 


EME 


H 


H 


C-H 
I 

h2c — o — ch2 — ch3 

DEE 


C; 


h3c 

l_l3C".  C-OH  +  H2C=CH2 

h3C  k  =  5.0  x  1 013  0.  63400.  S'1 

tert-  butanol  +  ethylene 

H 

I— 1,  I 

H2C=CH2  +  M"  C-OH 

H 

ethylene  +  methanol 

k  =  5.0  x  101?  0.  65200.  s1 


H2C=CH2  +  H3C-CH2— OH 


ethylene  +  ethanol 
Yasunaga  etal.  Comb.  Flame  (2011)  755/1032-1036. 

MP4/cc-pVTZ//MP2/cc-pVTZ  level  of  theory  with  zero  point  corrections 
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Ethylene  is  very  fast  to  ignite 
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3. 
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General  reaction  scheme 


C; 


High  temperature 


6® 


.0° 


o® 


_  Fuel - -  ^ 


/ 

\ 


„  ,  '(3-scission  ■  ^ Qo 

smaller  molecules  R 


H-atom  abstraction/ 

* 

5  "R^r 


olefins, 
cyclic  ethers,  <  i 
|3-scission 

Intermediate  temperature 


\7 


QOOH 


Low  temperature 

+  o2  chain 
branching 
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H-atom  abstraction  reactions 
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Acetone:  Low  p,  High  T 

1.25%  fuel,  22%  fuel  conversion,  T  =  1400  K,  p  =  1  atm 


CO  +  CH 


c2h6 


CH3CO  +  CH3 


/ 


53%  (f  =  0.5) 
64%  (f  =  1 .0) 
70%  (f  =  2.0) 


O  26%  (f  =  0.5) 


I 


c2h6 
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CH3COCH3 


1E12- 


1  El  1  - 


\ 


CHqCOCHo  +  HoO 


C; 


This  study 

Gierczak  et  al.  03  (PLIF) 
Wollenhaupt  et  al.  00  (PLIF) 
Yamada  et  al.  03  (PLIF) 
Vasudevan  et  al.  05 


”• - 1 - 

3 

1000  K/  T 


T 

4 


T 

5 
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(ts  JOU1  £UIO)  >1 


CHoOCH,  +  OH  =  CHoOCH,  +  H,0 


c; 


1E13- 


1  El  2  -i — r 
0.5 


Curran  et  al.  00 
Ogura  et  al.  07 
Bonard  et  al.  02 
Tranter/Walker  01 
Atkinson  97 
Arif  et  al.  97 
Mellouki  et  al.  95 
Nelson  et  al.  90 
Wallington  et  al.  88 
Tully/Droege  87 
Perry  et  al.  77 


1.0  1.5 


2.0  2.5  3.0 

1000  K/  T 
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k  (cm3  mol"1  s'1) 


CHoOCH,  +  OH  =  CHoOCH,  +  H,0 


c; 


1E13- 


1E12 


Cook  et  al.  09 
Curran  et  al.  00 
Ogura  et  al.  07 
Bonard  et  al.  02 
Tranter/Walker  01 
Atkinson  97 
Arif  et  al.  97 
Mellouki  et  al.  95 
Nelson  et  al.  90 
Wallington  et  al.  88 
Tully/Droege  87 
Perry  et  al.  77 


0.5  1.0  1.5 


2.0  2.5  3.0 

1000  K/  T 
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H-atom  Abstraction 


C 


>  n-Butanol  +  *OH/HO. 


C-W.  Zhou,  J.  M.  Simmie,  H.  J.  Curran.  Combust.  &  Flame,  2011  758  726-731. 

C-W.  Zhou,  J.  M.  Simmie,  H.  J.  Curran.  Int.  J.  Chem.  Kinet.,  2012  44  155-164. 

>  n-Butanol  +  CH3  (Imperial  College  London:  Prof.  Alex  Taylor ) 

D.  Katsikadakos,  C-W.  Zhou,  J.  M.  Simmie,  H.  J.  Curran,  P.A.  Hunt, 

Y.  Hardalupas,  A.M.K.P.  Taylor.  Proc.  Comb.  Inst.,  2012.  Paper  4D03  Thursday  2nd  August. 


>  Ethers  +  ‘OH 


EME 

C-W.  Zhou,  J.M.  Simmie,  H.J.  Curran 
Phys.  Chem.  Chem.  Phys.  2010  727221-7233. 
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H-atom  Abstraction 


>  Ketones  +  ‘OH 


C 


C-W.  Zhou,  J.  M.  Simmie,  H.  J.  Curran.  Phys.  Chem.  Chem.  Phys.  2011  13 11175-11192. 


>  /so-Butanol  decomposition  and  related  (CH3)2CH  +  CH2OH 
reaction  (Argonne  National  Laboratory:  Dr.  Stephen  Klippenstein) 


C-W.  Zhou,  S.  J.  Klippenstein,  J.  M.  Simmie,  H.  J.  Curran 
Proc.  Comb.  Inst.,  2012  in  press 
Paper :  4D05  Thursday  2nd  A  ugust  ([  0 
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log  (k/ cm3mol 


CH3OH  +  H02 


-  ch2oh  +  h2o2 


C 


2000 


77  K 


1000 


667 


500 


A 

n 

E 

Zhou 

6.24E-05 

4.89 

1.05E+04 

Klippenstein 

2.28E-05 

5.06 

1.02E+04 

Curran 

1.08E+04 

2.55 

1.05E+04 

Li 

3.98E+13 

0.00 

1.94E+04 

k  =  ATnexp(-E/RT) 

—  in  cm3,  mol,  s,  cal  units 


Klippenstein,  S  J.;  Harding,  L  B.;  Davis,  M  J.  et  ol.  Proc.  Comb.  Inst.  2011  33  351-358. 

Li,  J.;  Zhao,  Z.;  Kazakov,  A.;  Chaos,  M.;  Dryer,  F  L.  et  al.  Int.  J.  Chem.  Kinet.  2007  39  109-136. 
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Acetone:  a'  hydrogen  reactivity 


RCla 
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kl  cm3  mol'1 


Acetone:  a'  hydrogen  reactivity 


77  K 

2000  1000  667  500 


1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

0.5  1.0  1.5  2.0 

1000  K/  T 


>  Electron  delocalization  from  ti(CO)  to  a*(CHin)  and 
o*(CH0Ut)  is  different 


A£  =  2.0  kcal/mol 
AE  =  0.0  kcal/mol 
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k  /  cm3  mol'1 


a’  hydrogen  reactivity  in  DME 


c3 


77  K 

2000  1000  667  500 


AFct=  4.9  kcal/mol 
A£"ct=  1.7  kcal/mol 


>  Electron  delocalization  from  oxygen  lone  pairs  to 
a*(CH0Ut)  weaken  the  BDE  of  out-of-plane  CH  bond 
by  5.0  kcal/mol. 
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a'  Reactivity  Compared  to  Alkanes 


>  Ethers 


77  K 

2000  1000  667  500 


★  Oxygen  lone  pairs  accelerate  reactivity  of  a’  hydrogen  compared  to  alkane 

★  Growing  size  of  the  a-side  has  no  influence  on  reactivity  of  a’  hydrogen  atoms 


>  Ketones 


★  a  H-atom  is  less  reactive  than  primary  H-atom  in  alkanes 

★  Growing  size  of  the  a-side  will  accelerate  reactivity  of  a’  H-atoms 
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kl  cm  mol 


Comparative  Reactivity 


77  K 


2000  1000  667  500 


1000  K  IT  1000  Kl  T 
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General  reaction  scheme 


C 


High  temperature 


A© 


:o° 


O© 


Fuel 


H-atom  abstraction 


v 


Smaller  molecules  -  Mcission 


+  o. 


i>  ro2  ) 


olefins, 
cyclic  ethers,  <  ■ 

|B-scission 

Intermediate  temperature 


\7 


QOOH 


Low  temperature 

+  02 
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Low-temperature  chemistry 
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Alcohol  Oxidation 


Aldehyde  formation  from  a-radical  +  02 
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Alcohol  Oxidation 


Waddington  mechanism  (p-radical  +  02) 
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Comparison:  Alcohol/Alkane  Oxidation 


LLNL 


1000 


C 

bJO 


100 


"  10 

E 


0.1 


0.01 


♦  n-butanol  15  bar 
A  n-butane  18  bar 


♦♦ 


A 


"A  A" 


,'/  A  V 

^ ^ 


t£ 


>- 


'  / 


+ 


+ 


0.6 


0.8 


1  1.2 

1000/T [1/K] 


1.4 


1.6 
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Effect  of  chain  length 
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Alcohol  Oxidation 


C; 


Effect  of  chain  length  on  influence  of  functional  group 


n-Butanol  vs  n-Butane 


1000 


100 


”  10 

E 


0.1 


0.01 


\ 

- 

♦  n-butanol  15  bar 

A  n-butane  18  bar 

/ 

✓ 

: 

jf  / 

^  / 

1 

: 

: 

jf 

- 

V 

; 

- 1 - 1 - b- 

— ■ — i — ■ — i — ■ — i 

0.6  0.8  1  1.2  1.4  1.6 

1000/T [1/K] 


n-Pentanol  vs  n-Pentane 


1000 


■  n-pentanol  9  bar 
100  -H  •  n-pentane  8  -  11  bar 


„  1° 
if) 

E 


l 


o.i 


0.01 


■  • 


■  *■■■  • 

• 


■ . 


-h 


H- 


0.7  0.9 


1.1  1.3 

1000/T [1/K] 


1.5 
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Heufer  et  al.  Proc.  Comb.  Inst.  (2012)  in  press. 
Paper  4D06,  Thursday  2nd  August 

36 


NUI  Galway 
OE  Gaillimh 


Soybean  and  rapeseed  derived  biodiesels  have 

only  5  principal  components 


j 


R — ^  triglyceride  R 

o — \  o — / 

J  o  +  3CH30H 

methanol 


Fatty  acid  methyl  esters  (FAMEs): 

Methyl  Palmitate  (Cl 6:0) 


.OH 


.OH 


HO 


glycerol 


Methyl  Stearate  (Cl 8:0) 


Methyl  Oleate  (C18:l) 


Methyl  Linoleate  (Cl 8: 2) 


Methyl  Linolenate  (Cl 8: 3) 


C16:0  C18:0  C18:l  C18:2  C18:3 

Westbrook  et  al.  Proc.  Comb.  Inst.  (2012)  in  press. 
Paper  3D03  Wednesday  1st  August 
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Ignition  delay  -  ms 


Biodiesel  components  ignite  in  order  of  number  of| 
double  bonds 


C; 


LLNL 


Westbrook  et  al.  Proc.  Comb.  Inst.  (2012)  in  press. 
Paper  3D03  Wednesday  1st  August 
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C  =  C  double  bonds  reduce  low  T  reactivity 


s 

s 

a 

V  V 

a 

C-C- 

C- 

C-C 

-C 

s 

s 

a 

a 

s  s 
-C-C- 
s  s 


>  Inserting  one  C=C  double  bonds  changes  the  reactivity  of  4 
carbons  atoms  in  the  C  chain 


>  Allylic  C  -  H  bond  sites  are  weaker  than  most  others 

>  Therefore  they  are  preferentially  abstracted  by  radicals 

>  02  is  also  very  weakly  bound  at  allylic  sites  and  falls  off 
rapidly,  inhibiting  low  T  reactivity 


Westbrook  etal.  Proc.  Comb.  Inst.  (2012)  in  press. 
Paper  3D03  Wednesday  1st  August 


LLNL 
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Observed  same  effect  in  hydrocarbon  fuels:  hexenes 


Ignition  delay  times  in  a  rapid  compression 
machine  of  hexene  isomers 
(0.86-1.09  MPa,  0=1): 


C  =  C  -  C  -  C  -  C  -  C  1 -hexene 

C  -  C  =  C  -  C  -  C  -  C  2-hexene 

C  -  C  -  C  =  C  -  C  -  C  3-hexene 


R02  isomerization  initiates 
low  temperature  reactivity 

Moving  the  double  bond  towards  the  center 
of  the  molecule  “inhibits”  R02  kinetics 


Experimental  data:  Vanhove  et  al.  PCI  2005 
Simulations:  Mehl,  Vanhove,  Pitz,  Ranzi  Combustion 
and  Flame  155(2008)756—772. 
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Novel  fuels 
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Furans  vs  bio-ethanol 


■ 


>  Promising  Next  Generation  Biofuels 

>  Roman-Leshkov  et  al.,  Nature  (2007)  447: 
982-985.  (2nd  generation) 

>  Novel  renewable  production  process 

>  Biomass  (lignocellulosic)  feedstock  not 
destined  for  human/animal  consumption 

>  Highly  efficient 

>  Large  scale  and  low  cost 

>  Desirable  physicochemical  properties 

>  Higher  energy  density  (40%) 

>  Direct  combustion  in  unmodified  engine 

>  RON  =  119 

>  Lower  aqueous  solubility  and  less  volatile 


1  -i - ■ - 1 - ■ - 1 - ■ - 1 - ■ - 1 

0.8  1.0  1.2  1.4  1.6 


Laminar  Flame  Speed  Measurements  of  DMF,  Ethanol 
and  Gasoline,  Tian  et  al.,  Energy  Fuels,  2010 

Energy  Densities: 

>  Gasoline:  35  MJ/L 

>  2,5-DMF:  31  MJ/L 

>  Bio-Ethanol:  23  MJ/L 


Session  1 :  Monday  morning:  “Kinetics  of  Cyclic  Ethers” 
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2MF:  Unimolecular  Decomposition 


Energies  (0  K,  kJ  /  mol)  CBS-QB3,  CBS-APNO,  G3 

Somers  et  al.  Comb.  Inst.  (2012)  in  press.  PaperlDOl  Monday  30th 
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Conclusions 


>  General  chemical  reaction  schemes  of  HCs  can 
be  applied  to  oxygenated  fuels 

>  Details  of  oxygenated  fuel  combustion  are  quite 
different! 
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High-Pressure  Kinetic  Mechanisms  for 
Hydrogen  and  Hydrogen  Syngas 


1st  International  Workshop  on  Flame  Chemistry 

Warsaw,  Poland 
July  28,  2012 


Michael  P.  Burke 

Chemical  Sciences  and  Engineering  Division,  Argonne  National  Laboratory 

Frederick  L.  Dryer 

Department  of  Mechanical  and  Aerospace  Engineering,  Princeton  University 


Other  collaborators:  Yiguang  Ju,  Marcos  Chaos,  Jeffrey  Santner,  Francis  M.  Haas 

Stephen  Klippenstein,  Lawrence  Harding 


Motivation 


Growing  interest  in  computational 
engine  design/testing 
□  Fluid  mechanics  and  kinetics  sub-models 


H2  and  H2/CO 

□  Synthesis  gas  (H2/C0/H20/C02)  from 
coal/biomass  gasification 

□  Core  sub-model  for  all  fuels 


Advanced  engine  technologies 
High  P,  low  Tf 

□  Modeling  difficulties  for  flames 


(Shi  &  Reitz  2010) 


(Burke  et  at.  2011) 


1  ■  Y.  Shi,  R.D.  Reitz,  Fuel  89  (201 0)  341 6-3430. 

2.  M.P.  Burke,  M.  Chaos,  F.L.  Dryer,  Y.  Ju,  Combustion  and  Flame  157  (2010)  618-631. 

3.  M.P.  Burke,  F.L.  Dryer,  Y.  Ju,  Proceedings  of  the  Combustion  Institute  33  (201 1 )  905-912. 


2 


Mass  burning  rate  (g  cm 


Difficulty  in  predicting  high-pressure  flames 


Large  variations  among  models 

None  of  the  models  capture  pressure  dependence  across  all  conditions 


1  ■  M.P.  Burke,  M.  Chaos,  F.L.  Dryer,  Y.  Ju,  Combustion  and  Flame  157  (2010)  618-631 . 

2.  M.P.  Burke,  F.L.  Dryer,  Y.  Ju,  Proceedings  of  the  Combustion  Institute  33  (201 1 )  905-912. 
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H  mole  fraction 


What  controls 


Temoerature 


□ 


□ 


Sensitivity  Coefficient 


□ 


high-P/low-  Tj  flames? 

^^OH+O  (R1) 

H+02 

+M^  H02  (R2) 

increasing  P 

More  H02 

->  more  H02+radical  flux 

Flame  zone  shifts 
->  peak  sensitivity  at  higher  T’s 
->  collision  efficiencies  of  products 

More  R1/R2  competition 
->  amplified  sensitivity 

(Situation  similar  for  H2/CO) 


1  ■  M.P.  Burke,  M.  Chaos,  F.L.  Dryer,  Y.  Ju,  Combustion  and  Flame  157  (2010)  618-631 . 

2.  M.P.  Burke,  F.L.  Dryer,  Y.  Ju,  Proceedings  of  the  Combustion  Institute  33  (201 1 )  905-912. 
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Complexity  of  the  modeling  problem 


ho2+o=o2+oh 

H02+H02=H202+02 

o+h2=h+oh 

ho2+h=oh+oh 

ho2+oh=h2o+o2 

h2+oh=h2o+h 

h+o2=o+oh 


H2/02/He,  cp=0.30 
Tf  ~1400K 


H+02(+M)=H02(+M)| 


-2.5  -1.5  -0.5  0.5  1.5 

Sensitivity  Coefficient 


Uncertainty  in  all  reactions  of  10% 
-^burning  rate  uncertainty  of  30% 

Realistic  accuracy  improvements  for 
elementary  reactions  will  not  yield  typical 
expected  accuracies  for  global  behavior 

Optimization  against  global  targets 
necessary 


Functional  temperature  dependence  of 
0H+H02=H20+02  highly  disputed/ 
unknown 

Parameter  optimization  techniques  don’t 
work  if  the  functional  dependence  is  not 
known 


1  ■  M.P.  Burke,  F.L.  Dryer,  Y.  Ju,  Proceedings  of  the  Combustion  Institute  33  (201 1 )  905-912. 
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Complexity  of  the  modeling  problem 


ho2+o=o2+oh 

H02+H02=H202+02 

o+h2=h+oh 

ho2+h=oh+oh 

ho2+oh=h2o+o2 

h2+oh=h2o+h 

h+o2=o+oh 


H2/02/He,  cp=0.30 
Tf  ~1400K 


H+02(+M)=H02(+M)| 


-2.5  -1.5  -0.5  0.5 

Sensitivity  Coefficient 


1.5 


A  rigorous  modeling  solution  will  likely  require  both : 

□  Empirical  adjustments  to  rate  constants 

□  Improved  fundamental  understanding  of  select  processes 
Neither  alone  appears  sufficient  to  solve  the  problem. 


1  ■  M.P.  Burke,  F.L.  Dryer,  Y.  Ju,  Proceedings  of  the  Combustion  Institute  33  (201 1 )  905-912. 
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Updated  kinetic- transport  models 

■  H2:  Hong  et  al.  (201 1 )  and  Burke  et  al.  (2012)* 

□  H02  formation/consumption 

-  H+02(+M)  =  H02(+M) 

■  H02+radical  reactions 

□  H202  reactions 

□  ...  among  others 

■  CO:  Haas  et  al.  (2012) 

□  CO  +  OH  =  C02  +  H,  CO  +  H02  =  C02  +  OH 

□  HCO  chemistry 

*  Uncertainties  remained:  adjustments  of  rate  parameters  to 
improve  predictions 


1  ■  Z.  Hong,  D.F.  Davidson,  R.K.  Hanson,  Combust.  Flame  158  (2011)  633-644. 

2.  M.P.  Burke,  M.  Chaos,  Y.  Ju,  F.L.  Dryer,  S.J.  Klippenstein,  Int.  J.  Chem.  Kinet.  44  (2012)  444-474. 

3.  F.M.  Haas,  S.  Vranckx,  M.  Chaos,  R.X.  Fernandes,  F.L.  Dryer  (2012)  in  preparation. 
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Model  performance 


Time  (ms)  (Santer  et  al.  201 2, 5E01  on  Friday) 

■  Hong/Burke  perform  similarly  well  against  most  targets 

■  Largest  differences  in  flames 

□  Burke  et  al.  -  within  20%,  Hong  et  al.  -  within  40% 

■  Parameter  adjustments  not  unique  ->  uncertainties  remain! 

1  ■  Z.  Hong,  D.F.  Davidson,  R.K.  Hanson,  Combust.  Flame  158  (2011)  633-644. 

2.  M.P.  Burke,  M.  Chaos,  Y.  Ju,  F.L.  Dryer,  S.J.  Klippenstein,  Int.  J.  Chem.  Kinet.  44  (2012)  444-474. 

3.  F.M.  Haas,  S.  Vranckx,  M.  Chaos,  R.X.  Fernandes,  F.L.  Dryer  (2012)  in  preparation.  8 

4.  J.  Santner,  F.L.  Dryer,  Y.  Ju,  Proc.  Combust.  Inst.  (2012)  in  press,  oral  presentation  :  5E01  on  Friday. 


Uncertainties  remaining  in  2012  (for  flames) 

■  Parametric  uncertainties 

□  H02  +  X  reactions 

■  HOo  +  H  =  OH  +  OH 

Ho  +02 

h2o  +0 

■  ho2+oh  =h2o  +o2 

■  ho2  +  H02  =  H202  +  o2 

□  H  +  02  (+M)  =  H02  (+M) 

■  Pressure  dependence 

■  3rd  body  efficiencies  for  H20  and  COz 

□  CO  +  O  +  M  =  C02  +  M 

■  Model  assumptions 

□  Nonlinear  mixture  rules 

1  ■  M.P.  Burke,  M.  Chaos,  Y.  Ju,  F.L.  Dryer,  S.J.  Klippenstein,  Int.  J.  Chem.  Kinet.  44  (2012)  444-474. 

2.  F.M.  Haas,  S.  Vranckx,  M.  Chaos,  R.X.  Fernandes,  F.L.  Dryer  (2012)  in  preparation. 

3.  P.  Saxena,  F.A.  Williams,  7th  US  National  Combustion  Meeting,  Atlanta,  GA  ,  201 1 .  9 


Recall  the  complexity  of  the  modeling  problem 
and  uncertainties  in  0H+H02  =  H20+02 


A  rigorous  modeling  solution  will  likely  require  both : 

□  Empirical  adjustments  to  rate  constants 

□  Improved  fundamental  understanding  of  select  processes 
■  Neither  alone  appears  sufficient  to  solve  the  problem. 

1  ■  M.P.  Burke,  F.L.  Dryer,  Y.  Ju,  Proceedings  of  the  Combustion  Institute  33  (201 1 )  905-912. 
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Modeling  strategies 

Current  kinetic  models:  sets  of  rate  parameters 

□  Hierarchical,  comprehensive  modeling 

Westbrook  &  Dryer  (1984) 

□  Optimization  and  Uncertainty  Quantification 

Frenklach  (1984),  Frenklach,Wang,Rabinowitz  (1992):Solution-mapping  +  optimization  of  /4-factors 
Frenklach  et  al.  (2004),  Sheen  &  Wang  (2009):  Uncertainty  Quantification  of  /4-factors 

■  Turanyi  et  al.  (2012),  Sheen  et  al.  (2012):  Uncertainty  quantification  of  A-n-Ea 

□  Require  massive  amounts  of  data  to  constrain  full  7/P//lf  dependence  of  all  Ks 

■  Extrapolation  outside  the  dataset  very  challenging 

Direct  incorporation  of  theory  useful 

□  Replaces  fitting  formulas  with  physical  theories 

□  Common  for  extrapolation  of  data  for  a  single  reaction 

□  Imposes  constraints  spanning  all  T/P/M 

Multi-scale  models:  sets  of  molecular  parameters 

□  Optimal  use  of  information  from  ab  initio  calculations,  k  measurements, 
combustion  measurements 

□  Theory  fills  in  the  gaps  across  all  T/P/M 


1  ■  M.P.  Burke,  S.J.  Klippenstein,  L.B.  Harding,  Proceedings  of  the  Combustion  Institute  (2012)  in  press. 
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OH  +  H02  =  H20  +  O 


0.5  1  1.5  2  2.5  3  3.5  4 

1000/T  (K'1) 

Peeters  and  Mahnen  (1973) 

*  Hippier  etal.  (1995) 

▼ 

DeMore  (1979) 

Kappel  et  al.  (2002) 

Lii  et  al.  (1980) 

Srinivasan  et  al.  (2006) 

# 

Cox  et  al.  (1981) 

■  Hong  etal.  (2010) 

Kurylo  et  al.  (1 981 ) 

Keyser  (1 988) 

• 

Braun  et  al.  (1982) 

Sivaramakarishnan  et  al.  (2007) 

DeMore  (1982) 

“““■  Chaos  &  Dryer  (2008)  -  Hippier 

+ 

Goodings  &  Hayhurst  (1 988) 

. Chaos  &  Dryer  (2008)  -  Kappel 

□ 

Keyser  (1988) 

▲ 

Hippier  &Troe  (1992) 

(Rl) 

H202(+M)  =  OH+OH(+M) 

(R2) 

h2o2+oh  =  ho2+h2o 

(R3) 

HO2+HO2  —  H2O2+O 

2 

(R4) 

ho2+oh  =  h2o+o2 

(R5) 

OH+OH  =  0+H20 

12 


Multi- scale  informatics 


set  of  molecular  parameters  informed  by  data  across  all  scales 


I.  Molecular 
data 


II.  Rate  constant  III.  Combustion 

measurements  measurements 


Et  i/’s  v 

I—  ,  V  V  , 


imag’ 


TST, 

RRKM-ME,... 

- ► 


0-D  reactor, 
1-D  flame,... 


[OH]  vs.  t,s 


l/5 


^Mathematical  implementation  ^ 

□  Local  “surrogate  model” 

□  Least-squares  error  minimization 

^  □  Iterated  until  converged _  . 

(l)  +  (IV)^^=t  S„=s9 


(II) . 

(III) 


“  dX .. 


dF,  3lnfc 


dX, 


X  =  Optimization  parameters: 

Molecular  parameters  +  experimental  conditions 
Ef,  v’s,  vimag,  ...  +  T,  P  [MJ,... 

Fi(Xj)  =  Yu±Zi 


r 

ZJSij(Xj-Xj)  =  Yi±Zi^ 

j 

Xj,opt  and  Cx 

j 

J 

13 


Implementation  for  H202  system 


Optimization  variables 


10 


H202(+M)  =  OH+OH(+M) 
H202+0H  =  ho2+h2o 
ho2+ho2  =  h2o2+o2 
ho2+oh  =  h2o+o2 

OH+OH  =  0+H20 

Shock-heated  H202/H20/02/Ar 
Shock-heated  H20/02/Ar 
Shock-heated  H202/Ar 


A'(i)>n(1),  E(1) 

E  (2)  ,  V'all(2)  ,  v'tr(2)  ,  v'ss(2)  >  v'imag(2)  >  Ew(2)  ,  rf H2 02  ,  V  tS(2) 


(3) 


v'all(3)  >  v'tr(3)  ,  v'ss(3)  >  v’imag(3)  ,  Ew(3)  ,  Tj'jS(3) 


E  (4g)  ,  V  all(4)  ,  V  tr(4g),  V  ss(4g)  ,  V  imag(4g)  ,  Ew(4g j  ,  7J  TS(4g) 
E  (4e)>  v>TS(4e)  ,  v'tr(4e)  >  V'ss(4e)  >  rf TS(4e)  ,  f  ' VRCTST,c(4) 

E  (5g)  ,  V  all(5)  ,  V  tr(5g)  ,  V  ss(5g)  >  V  imag(5g)  ,  Ew(5g ) 

E  (5e)>  VTS(5e)  >  v'tr(5e)  ,  v'ss(5e) 

T'i,  P'i,  M'H202,o,i  ,  M'H20,o,i  ,  M'02,o,i 

T'i,  P'i,  M'H20,o,i  ,  o2,o, i  ,  M'H,o,i 

T'i,  P'i ,  M'H202,o,i  ,  O'  1,H202  ,  ^ 2,H202  ,  ^  1,H02  ,  ^ 2,H02 


O 

e  ,o” 


500 


Optimization  Targets 

I.  Molecular  data: 

ab  initio  calculations  (Klippenstein/Harding) 

II.  Rate  constant  measurements: 

see  paper 

III.  Combustion  measurements: 

OH(t),  H20(t)  Shock-heated  H^/Ar  (Hong  et  al.  2009,2010) 
OH(t)  Shock-heated  H20/02/Ar  (Hong  et  al.  2010) 

abs2i5nm(t)  Shock-heated  H^/Ar  (Kappel  et  al.  2002) 

1  ■  M.P.  Burke,  S.J.  Klippenstein,  L.B.  Harding,  Proceedings  of  the  Combustion  Institute  (2012)  in  press. 


\  OH  + 

ho2  =  h2o  +  o2 

q'  x  D  Qd  ^ 

• 

• 

p  X 

X 

DeMore  (1 980)  #* 

DeMore  (1 982)  9 

Keyser  (1988) 

<]  Lii  etal.  (1980) 

+  Cox  etal.  (1981) 

,•  Kappel  et  al.  (2002) 
Srinivasan  et  al.  (2006) 
Hong  etal.  (2010) 

Kurylo  etal.  (1981) 

Braun  etal.  (1982) 

1000  1500 

T(K) 


2000 


2500 
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Mole  fraction  (ppm)  Absorbance 


Different  interpretations 


for  OH+HO, 


Much  weaker  7-dependence 
(Secondary  reactions) 


Lower  magnitude 
(Arbitrary  H  atom  doping) 


1  ■  M.P.  Burke,  S.J.  Klippenstein,  L.B.  Harding,  Proceedings  of  the  Combustion  Institute  (2012)  in  press. 
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Mole  fraction  (ppm)  Absorbance 


Consistent  description  of  0H+H02  =  H20+02 


Single  description  consistent  with: 

1 .  Ab  initio  calculations 

2.  Low-  T  k  measurements 

3.  High-7" raw  global  data 

Milder  7-dependence 

□  Minimum  near  1200  K 


1  ■  M.P.  Burke,  S.J.  Klippenstein,  L.B.  Harding,  Proceedings  of  the  Combustion  Institute  (2012)  in  press. 
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Mole  fraction  (ppm)  Absorbance 


Consistent  description  of  0H+H02  =  H20+02 


Simultaneous  weighting  of  diverse  data  types 

□  Theory  guides  experimental  interpretations 

Raw  data  and  careful  documentation 
extremely  powerful 


1  ■  M.P.  Burke,  S.J.  Klippenstein,  L.B.  Harding,  Proceedings  of  the  Combustion  Institute  (2012)  in  press. 
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Mole  fraction  (ppm)  Absorbance 


Consistent  description  of  0H+H02  =  H20+02 


</> 


o 

E 


CO 


E 

o 


DeMore  (1980) 
DeMore  (1982) 
Keyser  (1 988) 

Lii  et  al.  (1980) 
Cox  et  al.  (1981) 
Kurylo  et  al.  (1981) 
Braun  et  al.  (1982) 


Hong  et  al.  (2012) 

>•  Kappel  et  al.  (2002) 
Srinivasan  et  al.  (2006) 
,■  Hong  etal.  (2010) 

■■  A  priori  model 
“  Constrained  model 


500  1000  1500 

T  (K) 


2000 


2500 


Z.  Hong,  K.-Y.  Lam,  R.  Sur,  S.  Wang,  D.F.  Davidson,  R.K.  Hanson 

“On  the  rate  constants  of  OH  +  H02  and  H02  +  H02:  A  comprehensive  study  of  H202  thermal 
decomposition  using  multi-species  laser  absorption.” 

Combustion  Symposium:  5D11 

M.P.  Burke,  S.J.  Klippenstein,  L.B.  Harding 

“A  quantitative  explanation  for  the  apparent  anomalous  temperature  dependence  of 
OH  +  H02  =  H20  +  02  through  multi-scale  modeling.” 


Combustion  Symposium:  4D09 


1  ■  M.P.  Burke,  S.J.  Klippenstein,  L.B.  Harding,  Proceedings  of  the  Combustion  Institute  (2012)  in  press. 
2.  Z.  Hong,  K.-Y.  Lam,  R.  Sur,  S.  Wang,  D.F.  Davidson,  R.K.  Hanson,  Proc  Combust  Inst  (2012)  in  press. 


Conclusions 

■  High-pressure  syngas  flames 

□  Emphasize  H02  pathways  +  collision  efficiencies  of  C02/H20 

□  Inherently  difficult  to  model 

■  Rigorous  modeling  solutions 

□  Empirical  adjustments  based  on  global  targets 

□  Improved  fundamental  characterization 

■  Uncertainties  remain  in  both  1)  model  parameters  and 
2)  model  assumptions 

■  Moving  forward 

□  Incorporation  of  theory  to  fill  in  the  gaps 

□  Raw  data  and  careful  documentation 

□  Characterization  of  non-idealities/uncertainties  in  experiments 
and  theory 
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Thank  you. 
Questions? 
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Determination  of  the  joint 
uncertainty  domain  of  rate  parameters: 

The  case  study  of  hydrogen  combustion 


T.  Turanyi,  T.  Nagy,  I.Gy.  Zsely,  T.  Varga, 
R.  Palvolgyi,  I.  Sedyo,  C.  Olm,  H.  J.  Curran 


Laboratory  for  Chemical  Kinetics 

Institute  of  Chemistry,  Eotvos  University  (ELTE), 

Budapest,  Hungary 


Uncertainty  of  rate  coefficients 

Definition  of  uncertainty  factor  t 

f(r ) = log  w(k°{T)/kmin  (r))=iog|0(/cniax(r)//c°  (r)) 

Calculation  of  the  variance  of  In  /(from  uncertainty  factor  f 
(assuming  3<7deviation  between  In  k?  and  In  /(max): 


a 


In  10 


The  databases  of  evaluated  combustion  kinetic  rate  parameters 
e.g.  Baulch  etal.,  2005 

-  recommended  temperature  independent  (constant)  uncertainty 
(usually  selected  from  values  /  =  0.1 ,  0.3,  0.5,  0.7,  and  1.0) 


OR 


-  uncertainty  is  defined  at  some  temperatures  (or  temperature  intervals) 


Uncertainty  parameter  fin  the  databases 


Features  of  uncertainty  factor  f\ 


>  widely  used  in  the  data  collections  /  data  evaluations 
(Warnatz,  Tsang,  Konnov,  Baulch,  NIST) 

>  available  for  most  reactions 

>  some  uncertainty  analysis  works  have  used  it 

>  important  input  information  in  the  mechanism  optimization  works 
of  Frenklach  et  al.  and  Wang  et  al. 

>  no  statistical  background,  subjective  assessment 

>  Usually  reflects  a  wider  uncertainty  than  the  quoted  experimental 
uncertainty  of  the  "best"  measurement 

>  temperature  dependence: 

•  quoted  to  be  temperature  independent  (50%) 

•  defined  at  several  (2-3)  temperatures  (50%) 


Arrhenius  plot  (In  k  vs.  1/7)  with  uncertain  parameter  A: 


uncertainty  of 
=>  the  rate  coefficient 
is  independent 
of  temperature 


Arrhenius  plot  (In  /cvs.  1/7)  with  uncertain  parameters  A  and  E : 


uncertainty  of 
=>  the  rate  coefficient 
depends 
on  temperature 


extended  Arrhenius  expression:  k(r)  =  at" exp(- e / rt) 


linearized  form: 


ln{A(r)}  =  ln{A}+  n  ■  ln{r}-  {E/r}- {t}' 

£  0 


-V - '  ^ 


k(6 


a  n  0 


Matrix-vector  form  of  the  linearized  Arrhenius  equation: 


p  :=[a  n  £ 


The  covariance  matrix  of  the  Arrhenius  parameters  and 
its  relation  to  the  uncertainty  of  the  rate  coefficient: 


2P  =  (p-p)(p-p)T= 


^an  ^  aP n 


r  <J  <7  r  <J  (7 

an  a  n  as  a  s 


< 7 ' 


r  <7  <7  r  <7  <7 

as  a  s  ns  n  s 


rne°n°e 

a2. 


=>  the  temperature  dependent  standard  deviation  of  k 
can  be  calculated  from  a  quadratic  form. 
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Relation  between  the  rate  coefficient  and  the 
covariance  matrix  of  the  Arrhenius  parameters 


For  the  3-parameter  Arrhenius  equation 

+\l<?a(Jn 1  n °-2rae(Ja(Jed''-'lrne(Jn 1  n  # 

For  the  original  (2-parameter)  Arrhenius  equation: 


If  the  variance  of  k  is  known  at  several  (at  least  6)  temperatures 
the  covariance  matrix  of  the  Arrhenius  parameters  can  be  calculated 

The  uncertainty  limits  of  the  Arrhenius  parameters  can  be 
determined. 

If  In  k  has  normal  distribution  o  the  Arrhenius  parameters 
have  joint  multivariate  normal  distribution. 


Reaction  O  +  N20  — >  NO  +  NO 
temperature  dependence  of  the  uncertainty  of  the 

rate  coefficient  (1 000K-4000K) 


k(T)  =  A  exp(-  E/RT) 
k($)  =  a  —  s-  0~l 
S' =  32,134 
e  = 13930 


<7a  =  0,355 
<Je  =  588 
ras  -  0,9045 


Reaction  O  +  N20  — ►  NO  +  NO 


■  :  -  c 

1.25E0 

6.28E-1 

■*  3,15E-1 
1.58E-1 
7,92E-2 
3,97E-2 
1.99E-2 
9.97E-3 
5,00E-3 


31,5 


a 

32,0 


8,0x1  O'4  ‘ 


15000 


31,0 

16000 
15000 
14000 


13000 


12000 
2,4x1  O'3 


2,0x1  O'3 


32,0 

a 


32,5 


33,0 


14000  ^ 

13000 


12000 


1 ,6x10 


J3 

(50 


4,0x1  O'4 


0,0  ' 
31,1 


Temperature  dependent 
1 D  normal 

probability  density  function  (pdf) 
of  the  rate  coefficient 
(1000  K  -  4000  K) 


Temperature  independent 
multivariate  joint  normal 
probability  density  function  (pdf) 
of  the  Arrhenius  parameters 
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Estimation  of  the  temperature  dependent  uncertainty 


For  most  of  the  reactions  the  data  evaluations  define 

-  temperature  independent  f  ( constant  /( 7)  function) 

OR 

-  a  verbally  defined  f(J)  function 

"f=  0.1  at  800  K  raising  to  0.2  at  2000  K  ” 


Reassessment  of  the  f(T)  functions  is  needed! 

Example: 

H  +  02  =  OH  +  O 


Reaction  H  +  02  — ►  O  +  OH 

All  experimental  and  theoretical  (TST)  results  from  the  NIST  database 


14 

13- 


O 

E 

E  11- 

o 


9- 


8 


0,0000  0,0002  0,0004  0,0006  0,0008  0,0010  0,0012  0,0014 


(1  /  7)  /  K 


-1 


3000  K 


800  K 


Ig(/r/  cm3mol'1s'1) 


Definition  of  kmin  and  kmax  ranges  of  uncertainty 


2500  K  ,  1 500  K 

(1/7 )  /  K 1 


Empirical  fi  7)  points  and  fitted  fCT)  function 


Nagy,  T.;  Turanyi,  T. 

Uncertainty  of  Arrhenius  parameters 
Int.  J.  Chem.  Kinet.,  43,  359-378(2011) 


Temperature  dependence  of  the  uncertainty  of  the 
rate  coefficients  of  the  most  important  H/O  reactions 


The  procedure  above  was  repetead  for  several  H/O  reactions 


Sedyo,  I;  Nagy,  T;  Zsely,  I;  Turanyi,  T 

Uncertainty  of  the  Arrhenius  parameters  of  important  elementary  reactions  of  the 
hydrogen-oxygen  system 
Proceedings  of  the  ECM,  2011 


Comprehensive  determination 
of  the  uncertainty  of  the  rate  parameters 

Required  features 

1 .  The  determination  should  take  into  account  all  available  information 
Direct  measurements:  determination  of  the  rate  coefficient  of 

a  single  reaction  step  at  given  temperature,  pressure  and  bath  gas 
Indirect  measurements:  ignition  delay  time,  laminar  flame  velocity, 
measured  concentration  profiles 

Theoretical  calculations  /  estimations:  uncertainty  is  critical 

2.  Statistics  based  (reproducible)  definition  of  the  uncertainty 

3.  Joint  uncertainty  for  all  rate  parameters: 

Arrhenius  parameters  A,  n,  E  for  each  reaction  step, 

3rd  body  efficiencies,  parameters  of  pressure  dependence 

4.  Extensions: 

thermodynamic  data:  in  accordance  with  the  ATcT  method 

other  parameters  of  simulations:  e.g.  transport  data 


Step  1 :  Collection  of  large  amount  of  indirect  measurements 

Aim:  collection  of  all  direct  and  indirect  measurement,  ever  measured, 

related  to  hydrogen  combustion 

Reality:  pre-selection  was  needed,  many  obsolate  measurements 
(before  about  1950)  were  not  considered 

>  ignition  delay  time  measurements  in  shock  tubes:  796  data  points  in  55  data 
series  from  17  publications 

>  ignition  delay  time  measurements  in  rapid  compression  machines  (RCMs): 

623  data  points  in  44  data  series  from  one  publication 

>  Laminar  flame  velocity  measurements:  622  data  points  in  64  data  series  from 
23  publications 

>  concentrations  measured  in  PSRs:  152  data  points  in  9  data  series  from  one 
publication 

>  concentration-distance  profiles  in  turbulent  flow  reactors:  214  data  points  in 
16  data  series  from  one  publication 

Total:  2407  data  points  in  188  data  series  from  43  publications. 

All  experimental  data  were  encoded  in  PrIMe  XML  data  format. 


PrIMe  database 


PrIMe  (Process  Informatics  Model)  database  of  Michael  Frenklach: 
www.primekinetics.org 

Open  database  that  contains  the  description  of 
several  hundred  indirect  combustion  measurements. 

Special  XML  data  format: 

defines  all  experimental  conditions  and  measurement  data. 


We  have  created  computer  codes  for 

-  encoding  literature  data  OR  other  experimental  data  in  PrIMe  XML 
format 

-  carrying  out  simulations  based  on  a  PrIMe  datafile 

(ignition  in  shock  tube,  laminar  stationary  flame,  tubular  reactor) 

-  calculating  the  deviations  between 
the  experimental  and  simulated  data 


PrIMe-format  encoding  of  a  shock  tube  experiment 


<commonProperties> 

<property  description=,,pressure  behind  reflected  shock  waves" 

label="P5"  name="pressure"  units="atm"> 
<value>64</value> 

<uncertainty  bound="plusminus"  kind="relative" 
transformation="1  ">0.01  </uncertainty> 

</property> 

<property  name="initial  composition"> 

<component> 

<speciesLink  preferredKey="H2"  primelD="s00009809"  /> 
<amount  units="mole  fraction">0.0033</amount> 
</component> 

<component> 

<speciesLink  preferredKey="02"  primelD="s00010295"  /> 
<amount  units="mole  fraction">0.00167</amount> 
</component> 

<component> 

<speciesLink  preferredKey="Ar"  primelD="s00000049"  /> 
<amount  units="mole  fraction">0.99503</amount> 
</component> 

</property> 

</commonProperties> 


Step  2:  Cross  testing  of  mechanisms  vs.  data 

Mechanisms: 


Group  1 :  hydrogen  combustion  mechanisms 

>6  Conaire  2004: 10  species,  21  reactions 

6  Conaire,  M.;  Curran,  H.  J.;  Simmie,  J.  M.;  Pitz,  W.  J.;  Westbrook,  C.  K.  Int.  J.  Chem.  Kinet.  2004,  36,  603-622. 

>Konnov  2008: 10  species,  33  reactions 

Konnov,  A.  A.  Combust.  Flame  2008, 152,  507-528. 

>Hong  2011:  10  species,  31  reactions 

Hong,  Z.;  Davidson,  D.  F.;  Hanson,  R.  K.  Combust.  Flame  2011, 158,  633-644. 

> Burke  2012:  1 1  species,  27  reactions 

Burke,  M.  P.;  Chaos,  M.;  Ju,  Y.;  Dryer,  F.  L.;  Klippenstein,  S.  J.  Int.  J.  Chem.  Kinet.  2012, 44,  444-474. 

Group  2:  originally  developed  for  wet  CO  combustion 

>Zsely  2005: 10  species,  32  reactions 

Zsely,  I.  G.;  Zador,  J.;  Turanyi,  T.  Proc.  Combust.  Instit.  2005,  30/1, 1273-1281. 

>Sun  2007: 1 1  species,  32  reactions 

H.  Sun,  S.l.  Yang,  G.  Jomaas,  C.K.  Law,  Proc.  Comb.  Inst.  2007,  31, 439-446. 

>CRECK2012: 11  species,  21  reactions 

Frassoldati,  A.;  Ranzi,  E.;  Faravelli,  T.  CRECK  modeling  Group  Hydrogen/CO  mechanism  version  1201. 

http://creckmodeling.chem.polimi.it/kinetic.html 

>Rasmussen  2008:  10  species,  30  reactions 

C.  L.  Rasmussen;  J.  Hansen;  P.  Marshall;  P.  Glarborg,  Int.  J.  Chem.  Kinet.  2008,  40,  454-480. 

>Keromnes  2012: 12  species,  33  reactions 

A.  Keromnes,  W.  K.  Metcalfe,  K.  A.  Heufer,  N.  Donohoe,  A.  K.  Das,  C.  J.  Sung,  J.  Herzler,  C.  Naumann,  P.  Griebel,  O.  Mathieu, 
M.  C.  Krejci,  E.  L.  Petersen,  W.  J.  Pitz,  H.  J.  Curran,  Combust.  Flame  (submitted) 

Group  3:  originally  developed  for  hydrocarbon  combustion 

>GRI  3.0: 10  species,  29  reactions 

Smith,  G  P;  Golden,  D  M;  Frenklach,  M;  Moriary,  N  W;  Eiteneer,  B;  Goldenberg,  M;  Bowman,  C  T;  Hanson,  R  K; 

Song,  S;  Gardiner,  W  C;  Lissianski,  V  V;  Qin,  Z  GRI-Mech  3.0  http://www.me.berkelev.edu/gri  mech/  23/1 1/2011/ 

>Li  2007: 1 1  species,  25  reactions 

Li,  J.;  Zhao,  Z.;  Kazakov,  A.;  Chaos,  M.;  Dryer,  F.  L.;  Scire,  J.  J.  J.  Int.  J.  Chem.  Kinet.  2007,  39, 109-136. 

>USC  2007: 10  species,  28  reactions 

H.  Wang;  X.  You;  A.  V.  Joshi;  S.  G.  Davis;  A.  Laskin;  F.  Egolfopoulos;  C.  K.  Law  USC  Mech  Version  II. 

High-Temperature  Combustion  Reaction  Model  of  H2/CO/C1-C4  Compounds.  http://ignis.usc.edu/USC  Mech  ll.htm 
(May  2007) 

>San  Diego  2011 : 1 1  species,  21  reactions 
San  Diego  Mechanism,  version  2011-11-22.  http://combustion.ucsd.edu 


Comparison  of  experimental  data  and  simulation  results 


Objective  function  Q : 


1^1 

2(p)=— y_L y 

n  tr  y  tr 


N-  f  tz  mod 


rr  (P)  -  ^ 


exp 


A2 


ij 


V 


a(^;xp) 


y 


y  ij 
lnyiJ 


if  G(yfp )  ~  constant 

if  C7(ln  yfp )  ~  constant 


Yij  measured/calculated  rate  coefficient  OR 

measured/calculated  ignition  time/flame  velocity 
in  data  point  j  of  data  series  / 

g  standard  deviation  of  the  measured  data 

N,  number  of  data  points  in  data  series  / 

N  number  of  data  series  (different  experiments) 
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Performance  of  the  mechanisms  vs.  years 


Results  of  cross  testing 


average  of  error  function  values 

Ignition 

delays 

PSR 

flame 

Overall  performance 

Mechanism 

Asaba  etal.  (1965) 
only 

48.13 

26.34 

32.67 

43.34 

GRI  3.0  1999 

184.31 

13.89 

49.40 

30.80 

20.42 

6  Conaire  2004 

122.88 

21.87 

83.19 

54.34 

33.64 

Zsely  2005 

143.30 

14.00 

54.57 

32.32 

21.31 

Li  2007 

154.69 

17.87 

65.83 

29.27 

24.94 

Sun  2007 

159.96 

16.33 

51.98 

32.06 

22.70 

USC  2007  Version  II 

161.65 

15.12 

82.59 

31.52 

25.12 

Konnov  2008 

153.68 

31.55 

82.11 

36.14 

37.87 

Rasmussen  2008 

152.07 

16.80 

54.48 

29.78 

22.97 

Hong  2011 

165.11 

20.77 

55.60 

33.65 

26.61 

San  Diego  2011 

181.05 

18.69 

52.15 

31.66 

24.39 

Burke  2012 

171.97 

11.88 

57.08 

48.03 

22.43 

CRECK  2012 

157.02 

14.39 

65.39 

30.74 

22.56 

Keromnes  2012 

163.74 

634 

136 

200 

970 

No.  data  points 

37 

53 

8 

19 

80 

No.  datasets 

3 

Some  indirect  experimental  data  series  are  presumably  wrong  and  these 
were  excluded  from  the  optimization  and  mechanism  testing. 


example:  T.  Asaba,  W.C.  Gardiner,  R.F.  Stubbelman, 
Proc.  Combust.  Inst.,  10,  295-302(1965) 


Step  3:  Optimization  of  the  rate  parameters 


1.  Direct  and  indirect  measurements  are  taken  into  account 
simultaneously. 

2.  The  uncertainty  domain  of  all  the  Arrhenius  parameters  is  determined 


3.  For  each  important  reaction  step,  the  Arrhenius  parameters  A,  n,  E  and 
maybe  other  rate  parameters  are  determined 
within  their  domain  of  uncertainty. 


Turanyi  T,  Nagy  T,  Zsely  IGy,  Cserhati  M,  Varga  T,  Szabo  B, 

Sedyo  I,  Kiss  P,  Zempleni  A,  Curran  H  J 

Determination  of  rate  parameters  based  on  both  direct  and  indirect  measurements. 
Int.  J.  Chem.  Kinet.  44,  284-302  (2012) 


Optimization  and  uncertainty  estimation  of  parameters 


Optimization  =  minimization  of  this  error  function: 


yk  measured/calculated  rate  coefficient  OR 
measured/simulated  indirect  measurement 


Calculation  of  the  covariance  matrix  of  the  estimated  parameters: 


ijoW£yX)~1JoWEy‘](£T+Ej[(jIW£TlJ0)rljTW£ 


-1 

Y 
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Transformation  of  the  covariance  matrix  of  parameters 
to  ..traditional”  gas  kinetics  measures  of  uncertainty 

covariance  matrix  of  Arrhenius  parameters  covariance  matrix  of  In  k 


cov(k,  (T),  Kj  (T))  -  (K,  (T)  -  Kj  (T))(Kj  (T)  -  Kj  (T))  -  ©T  (Pi  -  Pi  )(Pj  -  Pj  )T0  =  ©%  _Pj  © 


0T  :=  [l  ln6  -0”1] 


1  1  A 

o(ln  k)  =  In  10  o(log10  k)  =  — —  f  (t) 

mf 


covariance  matrix  of  In  /(—» correlation  /"between  the  In  /rvalues 

~  cov(Ki  (t),  Kj  (t)) 

r-.(T)=  O  (T) G  (T) 

1  J 


Strategy  for  optimization 


STEP1 


STEP2 


STEP3 


experiment 

number  of  datapoints 

H+02=0+0H 

LPH+02(+M)=H02(+M) 

0+H2=H+OH 

0H+H2=H+H20 

H202+H=H2+H02 

H02+H=20H 

H02+H=H2+02 

H02+0H<=>H20+02 

Herzler  et  al.  (2009a) 

9 

O 

0 

Fujimoto  and  Suzuki  (1967) 

9 

0 

0 

Zhang  et.  al.  (201 2a) 

7 

0 

0 

Naumann  et.  al.  (2011 a) 

19 

0 

0 

Naumann  et.  al.  (201 1 b) 

26 

0 

0 

Petersen  et  al.  (2003a) 

9 

0 

0 

Cheng  and  Oppenheim  (1984a) 

58 

0 

0 

Petersen  et  al.  (2003b) 

24 

0 

0 

Petersen  et  al.  (2003c) 

4 

0 

0 

Petersen  et  al.  (1996a) 

16 

0 

0 

0 

Petersen  et  al.  (1996b) 

6 

0 

0 

0 

Slack  (1977) 

12 

0 

0 

0 

Bhaskaran  et  al.  (1973) 

14 

0 

0 

0 

Wang  et  al.  (2003a) 

12 

0 

0 

0 

Naumann  et.  al.  (2011°) 

13 

0 

0 

0 

Chaumeix  et  al.  (2007a) 

5 

0 

0 

0 

Chaumeix  et  al.  (2007b) 

7 

0 

0 

0 

Chaumeix  et  al.  (2007c) 

5 

0 

0 

0 

Cohen  et.  al.  (1967) 

21 

0 

0 

0 

Cheng  and  Oppenheim  (1984b) 

56 

0 

0 

0 

0 

Naumann  et.  al.  (2011 d) 

19 

0 

0 

0 

0 

Naumann  et.  al.  (2011 e) 

19 

0 

0 

0 

0 

Petersen  et  al.  (2003a) 

9 

0 

0 

Cheng  and  Oppenheim  (1984a) 

58 

0 

0 

Petersen  et  al.  (2003b) 

24 

0 

0 

Petersen  et  al.  (2003c) 

4 

0 

0 

Petersen  et  al.  (1996a) 

16 

0 

0 

0 

Petersen  et  al.  (1996b) 

6 

0 

0 

0 

Slack  (1977) 

12 

0 

0 

0 

Bhaskaran  et  al.  (1973) 

14 

0 

0 

0 

Wang  et  al.  (2003a) 

12 

0 

0 

0 

Naumann  et.  al.  (2011°) 

13 

0 

0 

0 

Chaumeix  et  al.  (2007a) 

5 

0 

0 

0 

Chaumeix  et  al.  (2007b) 

7 

0 

0 

0 

Chaumeix  et  al.  (2007c) 

5 

0 

0 

0 

Cohen  et.  al.  (1967) 

21 

0 

0 

0 

Cheng  and  Oppenheim  (1984b) 

56 

0 

0 

0 

0 

Naumann  et.  al.  (2011 d) 

19 

0 

0 

0 

0 

Naumann  et.  al.  (2011 e) 

19 

0 

0 

0 

0 

Zhang  et.  al.  (201 2b) 

10 

0 

0 

0 

Wang  et  al.  (2003b) 

10 

0 

0 

0 

Wang  et  al.  (2003c) 

21 

0 

0 

0 

Naumann  et.  al.  (2011 f) 

9 

0 

0 

0 

0 

Wang  et  al.  (2003d) 

12 

0 

0 

0 

0 

Naumann  et.  al.  (2011a) 

10 

0 

0 

0 

0 

Petersen  et  al.  (1996c) 

8 

0 

0 

0 

0 

Herzler  et  al.  (2009b) 

9 

0 

0 

0 

0 

Petersen  et  al.  (1996d) 

3 

0 

0 

0 

0 

0 

Zhang  et.  al.  (2012°) 

8 

0 

0 

0 

0 

Herzler  et  al.  (2009c) 

12 

0 

0 

0 

0 

Naumann  et.  al.  (201 1 h) 

16 

0 

0 

0 

0 

Petersen  et  al.  (1996e) 

14 

0 

0 

0 

Petersen  et  al.  (1996f) 

7 

0 

0 

0 

0 

Schott  and  Kinsey  (1958) 

17 

0 

0 

0 

0 

Petersen  et  al.  (1996a) 

17 

0 

0 

0 

0 

Naumann  et.  al.  (2011') 

18 

0 

0 

0 

0 

0 

Naumann  et.  al.  (20110 

13 

0 

0 

0 

0 

0 

0 

0 

Results  of  optimization 

A  new  set  of  optimized  rate  parameters 

(in  this  case:  26  optimized  parameters  (23  Arrhenius  parameters  of  8  reactions 

3  third  body  efficiencies) 

Covariance  matrix  of  all  optimized  parameters 
CONVERTED  TO 

f(T)  uncertainty  function  for  each  rate  coefficient 
r(T)  correlation  function  of  pairs  of  rate  coefficients 


2400  2000  1400  1200  1000 


literature  k  value 
literature  /(uncertainty 

optimized  /(value 
optimized  /(uncertainty 


reaction  0H+H2=H20+H: 


Step  4:  Testing  the  optimized  mechanism  against 

other  mechanisms 

The  optimization  code  created  figures  on 

the  agreement  of  the  simulation  results  and  the  experimental  points. 


c n 

O) 

10 


■  exp. 

- GRI  3.0  1999 

- Conaire  2004 

- Zsely  2005 

- Li  2007 

- Sun  2007 

- USC  2007  II 

- Konnov  2008 

- Rasmussen  2008 

- Hong  201 1 

- San  Diego  201 1 

- Burke  2012 

- CRECK  2012 

- Keromnes  2012 

optimized  mechanism 


0.90 


0.92 


0.94  0.96 

1000  K/  T 


0.98 


1.00 


experiments:  Pang,  G.A.,  Davison,  D.F.,  Hanson,  R.K., 
Proc.  Combust.  Inst.,  32,  181-188(2009) 

Fig.  3.,  full  square;  4%  H2  /  2%  02  /  Ar,  Pf=0  =  3.5  atm 


Testing  the  optimized  mechanism  against  other  mechanisms 


mechanisms 

set  #1 

GRI  3.0  1999 

43.63 

6  Conaire  2004 

11.01 

Zsely  2005 

16.64 

Li  2007 

11.47 

Sun  2007 

15.72 

use  2007  II 

14.31 

Konnov  2008 

12.42 

Rasmussen  2008 

28.58 

Hong  2011 

14.44 

San  Diego  2011 

18.18 

Burke  2012 

16.49 

CRECK  2012 

9.28 

Keromnes  2012 

11.94 

optimized 

mechanism 

8.25 

error  function  values  (smaller  is  better) 


set  #1 :  Ignition  time  experiments  used  at  the  optimization 


Testing  the  optimized  mechanism  against  other  mechanisms 


mechanisms 

set  #1 

set#2 

GRI  3.0  1999 

43.63 

48.13 

6  Conaire  2004 

11.01 

13.89 

Zsely  2005 

16.64 

21.87 

Li  2007 

11.47 

14.00 

Sun  2007 

15.72 

17.87 

use  2007  II 

14.31 

16.33 

Konnov  2008 

12.42 

15.12 

Rasmussen  2008 

28.58 

31.55 

Hong  2011 

14.44 

16.80 

San  Diego  2011 

18.18 

20.77 

Burke  2012 

16.49 

18.69 

CRECK  2012 

9.28 

11.88 

Keromnes  2012 

11.94 

14.39 

optimized 

mechanism 

8.25 

11.59 

set  #1 :  Ignition  time  experiments  used  at  the  optimization 
set  #2:  All  available  ignition  time  experiments 


Testing  the  optimized  mechanism  against  other  mechanisms 


mechanisms 

set  #1 

set#2 

set#3 

GRI  3.0  1999 

43.63 

48.13 

43.34 

6  Conaire  2004 

11.01 

13.89 

20.42 

Zsely  2005 

16.64 

21.87 

33.64 

Li  2007 

11.47 

14.00 

21.31 

Sun  2007 

15.72 

17.87 

24.94 

use  2007  II 

14.31 

16.33 

22.70 

Konnov  2008 

12.42 

15.12 

25.12 

Rasmussen  2008 

28.58 

31.55 

37.87 

Hong  2011 

14.44 

16.80 

22.97 

San  Diego  2011 

18.18 

20.77 

26.61 

Burke  2012 

16.49 

18.69 

24.39 

CRECK  2012 

9.28 

11.88 

22.43 

Keromnes  2012 

11.94 

14.39 

22.56 

optimized 

mechanism 

8.25 

11.59 

19.35 

set  #1 :  Ignition  time  experiments  used  at  the  optimization 

set  #2:  All  available  ignition  time  experiments 

set  #3:  All  types  of  data  (ignition  time,  PSR,  laminar  flame  velocity) 


Testing  the  optimized  mechanism  against  other  mechanisms 


mechanisms 

set  #1 

set#2 

set#3 

GRI  3.0  1999 

43.63 

48.13 

43.34 

6  Conaire  2004 

11.01 

13.89 

20.42 

Zsely  2005 

16.64 

21.87 

33.64 

Li  2007 

11.47 

14.00 

21.31 

Sun  2007 

15.72 

17.87 

24.94 

use  2007  II 

14.31 

16.33 

22.70 

Konnov  2008 

12.42 

15.12 

25.12 

Rasmussen  2008 

28.58 

31.55 

37.87 

Hong  2011 

14.44 

16.80 

22.97 

San  Diego  2011 

18.18 

20.77 

26.61 

Burke  2012 

16.49 

18.69 

24.39 

CRECK  2012 

9.28 

11.88 

22.43 

Keromnes  2012 

11.94 

14.39 

22.56 

optimized 

mechanism 

8.25 

11.59 

19.35 

set  #1 :  Ignition  time  experiments  used  at  the  optimization 

set  #2:  All  available  ignition  time  experiments 

set  #3:  All  types  of  data  (ignition  time,  PSR,  laminar  flame  velocity) 


The  future 


Detailed  hydrogen  combustion  mechanisms 


Based  on  the  current  trend, 

the  absolutely  accurate  hydrogen  combustion  mechanism 

will  be  available  on  7th  March,  2032  © 


Summary 

1  The  f{T)  uncertainty  parameters  available  from  the  data  collections 
are  realistic,  but  have  several  shortcomings,  like  it 

refers  to  the  uncertainty  of  k  only. 

Needed:  a  joint  uncertainty  of  the  rate  parameters. 

2  Assessment  of  the  joint  uncertainty  domain  of  the 
Arrhenius  parameters  of  a  selected  reaction  step 

based  on  all  direct  measurements  and  TST  calculations. 

•  upper  estimation  of  the  uncertainty 

•  no  information  about  the  correlation  of  two  /rvalues 

3  Comprehensive  evaluation  of  the  uncertainty  of  all  rate  parameters: 

For  a  given  chemical  kinetic  system, 

all  direct,  indirect  and  theoretical  results  are  collected 

assessment  of  the  uncertainty  of  input  data 

=>  new,  physically  meaningful  rate  parameters 
=>  information  about  the  joint  uncertainty  of  rate  parameters 

=>  better  simulation  results 

=>  quantification  of  the  uncertainty  of  the  simulation  results 
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The  role  of  uncertainty  quantification 
within  kinetic  model  improvement  in 
combustion 


Alison  S.  Tomlin 


Development  of  complex  kinetic 
mechanisms 


UNIVERSITY  OF  LEEDS 


In  combustion  complex  chemical  mechanisms  built  by: 

>  proposing  sets  of  rules  for  species  interactions: 

mechanism  protocol. 

>  developing  effective  parameterisations  for  the  kinetics 
described  within  the  mechanism. 

•  Ability  to  specify  protocols  well  developed  e.g.  hydrocarbon 
oxidation. 

•  Large  comprehensive  mechanisms  e.g.  biodiesel  surrogate 
methyl  decanoate:  3012  sp.,  8820  reacs  (Herbinet  et  al.,  2008) 

•  Many  rate  coefficients  have  to  initially  be  estimated  using  rules 
related  to  chemical  structure. 

•  Does  this  lead  to  a  robust  mechanism?  How  can  the  model 
robustness  be  improved? 


g2(a,e) 

CO  -> 


Improving  model  robustness 


Joint  pdf  of 
inputs 


UNIVERSITY  OF  LEEDS 


Probability  distribution 
of  predicted  Y, 


Model  response 

—T 


We  would  like  to  narrow  the  output  distribution  as  far 
as  possible  i.e.  to  reduce  predictive  error  bars. 


standard  deviation  /  ms 


Sensitivity  and  uncertainty 
analysis 


*  Uncertainty  analysis  (UA)  puts 
error  bars  on  predictions. 


•  Overlap  (or  not)  with 
experimental  values  tells  us 
whether  our  model  may  be 
structurally  robust. 


absolute  mean  perturbation  /ms 


UNIVERSITY  OF  LEEDS 


Sensitivity  analysis  (SA)  tries  to 
determine  how  much  each  input 
parameter  contributes  to  the  output 
uncertainty  (usually  variance). 


rrr 


Global  sensitivity/uncertainty  methods  university  of  leeds 


Global  -  attempts  to  cover  whole 
input  space  using  a  sampling 
method. 


Why  global  methods? 


rrr 


UNIVERSITY  OF  LEEDS 


Sampling  methods  can  be  computationally  costly  so  why 
bother. 

•  To  simulate  output  distributions  for  available  targets. 

•  To  calculate  the  contribution  of  important  parameters  to 
these  distributions  however  large  their  uncertainty  ranges. 

•  To  see  how  each  target  helps  to  constrain  a  parameter 

>  Could  be  important  for  the  use  of  optimisation  approaches  where 
parameter  isolation  helps  to  narrow  constraints. 

•  To  explore  parameter  interactions. 


Sampling  approaches 


rrr 


UNIVERSITY  OF  LEEDS 


•  Need  efficient  sampling 

approach  where  the  output 
moments  converge  quickly  with 
sample  size.  1 

•  Quasi-random  approaches  ° 
usually  work  best. 

•  Sample  from  input 
distributions  until  predicted  a| 
output  mean  and  variance 
converge. 

•  Then  need  a  method  to  i 
investigate  input-output  | 
relationships. 


Random  sampling  LH  sampling 


Halton  sequence  Sobol  sequence 


High  Dimensional  Model 
Representations  (HDMR) 


•  Input-output  relationship  is  expressed  as  a  finite  hierarchical  function 
expansion: 


n 


•  Model  replacement  built  using  quasi  random  sample  and  approximation  of 
component  functions  by  orthonormal  polynomials. 

•  Equivalent  to  ANOVA  distribution  so  each  term  represents  contribution  to 
overall  output  variance. 

•  1st  &  2nd  order  sensitivity  indices  easily  calculated  from  polynomial 
coefficients. 


Monte  Carlo  (MC)  simulations 


•  Component  functions  reveal  independent  mean  effect  of  each 
parameter  and  highlight  any  non-linear  responses. 

•  Local  S’;  at  the  nominal  value  may  be  quite  different  from  the 
overall  effect  across  the  feasible  parameter  range. 

•  Scatter  reveals  influence  of  other  parameters. 


Example  from 
flame  calculation: 
NOx  prediction. 


A-factor  (cm3  mol  1  s  ^  xio 


1st  order  sensitivities  and  B1 

ranking  +  component  functions  university  of  leeds 


Results 


1  st-order 


2nd-order 


Plots 


1  st-order 


2nd-order 


Output: 


1 


Exit 


Ranking 

1. 

2. 

3. 

4, 

5. 

Input 

47 

102 

138 

8 

142 

Si 

0.2619 

0.1 328 

0.1047 

0.0506 

0.0327 

Sum  Si  =  0.81 88 


Polynomial  Order : 
Sobol  Index : 
rA2  Pearson : 
rA2  Spearman : 


0.2619(1 .) 
0.2488  (1 .) 
0.2313(1 .) 


Tells  us  fraction  of  output 
variance  due  to  each  of 
the  inputs. 


8  0.02 


0.01 


1 02 


Polynomial  Order : 
Sobol  Index : 
rA2  Pearson : 
rA2  Spearman : 


0.1328  (2.) 
0.1234  (2.) 
0.1588  (2.) 


Input: 


102 


Scatter 


HDMR 


Both 


2nd  order  component  functions 
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Input  1 : 
4 


Input  2: 
5 


HDMR 


Polynomial  Order :  2 

Sobol  Index :  0.0124(3.) 


Polynomial  Order :  2 

Sobol  Index:  0.0145(2.) 


Typical  practice 
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•  Constructed  mechanism  incorporated  within  model  which 
simulates  chosen  experimental  targets: 

>  flame  speeds,  ignition  delays,  JSR  concentration  profiles 

•  Agreement  with  targets  assessed  and  sensitivity  analysis 
performed. 


Updates  made  to  key 
parameters  e.g. 
through  theory  or 
direct  experiments 


Mechanism  optimised 
using  available  targets 
from  indirect 
experiments 


rrr 
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Example  1 
Methanol  oxidation 


Comparison  of  simulated  ign.  delays  (Li  et  B 

al.  2007)  with  data  (Kumar  &  Sung,  201 1 )  university  of  leeds 


(a)  Methanol/02/Ar,  $=1.0,  bar 


Constraints  from  the 
experiments? 


8? 


L8 

L6 

L4 

L2 

I 

10.8  10.9  11  11.1  11.2 

logl0[/t  (cm3  moH  s-1}] 


Base  (A,  t) 
Experiment  (t) 


,t)  ' 


Line:  linear  le.r-u-  -qua re  Lit'1 

j _ i _ i _ i _ i _ L 


0.30 

0.25 


Impact  of  modifying  rate 
to  Tsang  value 


>. 


2  0.20 


o 

•0 

>, 


£  0.15 


-C 

-S 

o 


0.10 


0.050 


0.0 
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•  Mechanism  -  Li  et  al.  (2007);  18  species,  93  reactions. 

•  Target  output  -  ignition  delay  time  (t)  for  stoichiometric  mixtures  of 
methanol  and  oxygen  over  a  range  of  temperatures  and  pressures. 

•  Enthalpies  of  formation  and  ^-factors  varied  over  random  sample. 

•  Using  initial  ranges  one  reaction  dominates  (up  to  90%  of  total  output 
variance)  with  strong  nonlinearities. 


0.0006-1 


CH3OH+HO2 
Stage  1 


Low  scatter  indicates  low 
influence  of  all  other  parameters. 


0.0004  - 


0.0000 


0.0 


0.5 


1.0 


X 


(T,P,<p)=(1150K,5bar,1) 


H 

Sensitivity  mapping  over  PIT  university  of  leeds 


•  Can  also  test  how  sensitivities  change  outside  of  conditions  where 
experimental  data  exists  -  although  uncertainties  in  activation  energies 
were  not  included  in  this  study. 

•  No  dramatic  changes  in  this  case. 


r  0.9 


-  0.8 


-  0.7 


.  1100 
^  1200 
" Ofa  1300 


0  1 
0.0 


0 


k(T)  (cm3/mol/s) 


Model  Updates 
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•  CH30H+H02  a  key  reaction  with  initial  uncertainty  (1  O')  of  5. 

•  Transition  state  theory  (TST)  and  variational  TST  (VTST)  then 
used  to  re-estimate  its  rate  coefficient  (Klippenstein,201 1). 

•  Lower  than  value  in  Li  at  low  T  (8  times  i.e.  larger  than  5!)  with 
an  estimated  uncertainty  of  2. 


•  Close  agreement 
between  new  value  and 
that  of  Tsang  at  low 
temperatures. 

•  Closer  to  Li  value  at 
high  temperatures. 


Scaled  Probability 


H 

Results  of  model  updates  university  of  leeds 


t/<t> 


(T,P,<p)=(1150K,5bar,1) 


Further  impacts? 
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•  Modified  rates  from  both  studies  are  similar  at  low- 
intermediate  temperatures  -  well  constrained? 


•  Low  sensitivity  of  high  T  shock  tube  simulations  (Bowman 
1 975)  and  predicted  flame  speeds. 


•  High  sensitivity  to 
flow  reactor  data  at 
970  K  and  2.5  atm. 


CH30H+H02  =  CH20H+H202  _ 
H202(+M)  =  20H(+M)  — 

ch3oh+oh  =  ch3o+h2o 
ch3oh+oh  =  ch2oh+h2o 

H+02  =  O+OH 
HCO+M  =  H+CO+M 
OH+CH3(+M)  =  CH3OH(+M) 

ch2oh+m  =  ch2o+h+m 

CO+OH  =  C02+H 

ch3oh+h  =  ch2oh+h2 
hco+o2  =  co+ho2 
ch2o+oh  =  hco+h2o 
ho2+oh  =  h2o+o2 

2H02  =  H202+02 

1  1  I  1 

-0.6 


c 


I  |  Flame  Speed 
I  |  Ignition  Delay 
■■  VPFFT1000 


“i — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1— 

-0.4  -0.2  0.0  0.2  0.4  0.6 


b) 


Sensitivity  Coefficient 


How  to  combine  knowledge? 
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•  The  reaction  exhibits  sensitivities  for  several  reaction 
systems. 

•  Should  start  to  become  parameter  focussed  rather  than 
mechanism/target  focussed? 

•  Possibility  to  incorporate  all  sets  of  data  to  which  a 
parameter  is  sensitive  within  a  Bayesian  framework. 

•  Does  this  get  too  messy  for  complex  fuels? 


Uncertainty  factor, 


Possible  for  Ho/O. 
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0.6 
0.5  b 
0.4 

0.3 

0.2 

0.1  - 

°0 


♦ — », 

>  V 


500 


Baulch  et  al. 


**.,  A  priori  model 


H02+0H=H20+02  Burke  et 
al.  2012 

Data  from  TST,  discharge 
flow,  shock  tubes  used. 


Constrained  model 


1000  1500  2000  2500 

T  (K) 


H+02=0H+0, 
Turanyi  et  al.  2012 

Data  from  ignition 
delays,  shock 
tubes 


Larger  fuels? 
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•  Target  predictions  for  more  complex  mechanisms  are 
unlikely  to  be  driven  by  only  1  or  2  parameters. 

•  What  type  of  approach  should  be  adopted? 

•  A  hierarchical  one  where  the  simpler  and  better  constrained 
systems  are  tackled  first? 

>  Turanyi’s  approach  only  selects  data  which  is  strongly  sensitive  to 
parameters  being  constrained  i.e.  reasonably  well  isolates 
parameters. 

•  The  resulting  more  tightly  constrained  parameters  could 
then  be  included  in  studies  of  more  complex  fuels  -  helping 
to  isolate  the  other  reactions. 


rr 


UNIVERSITY  OF  LEEDS 


Example  2:  Butane  oxidation  in 
a  jet  stirred  reactor 


Performance  of  EXGAS 
mech. 


Isothermal  jet-stirred 
reactor 

1  atmosphere 

Residence  time  :  6  s 

Equivalence  ratio  of  1 

4%  butane  as  inlet 
mole  fraction. 


IQ- 


550 
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600  650  700  750  800 

T/K 


Methodology 
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•  EXGAS  mechanism  -  1304  uncertain  A-factors  for  forward 
reactions  studied. 

•  Very  large  dimensional  input  parameter  space. 

•  Screen  out  unimportant  parameters  before  full  global 
analysis. 

•  Linear  sensitivities  at  several  temperatures  used  for 
screening. 

•  Gives  better  coverage  for  remaining  parameters  allowing 
smaller  sample  size. 


Global  Analysis 
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•  31  /4-factors  selected  for  global  runs  based  on  the  linear 
screening. 

•  Input  uncertainty  factors  assigned  from  appropriate 
evaluation  studies  where  available  -  C0-C2  reactions  base 
(Baulch  et  al.  1 994) . 

•  For  calculated  parameters  within  the  primary  and 
secondary  mechanisms,  10f  ~  3. 

•  n-butane  +  OH  1 0f  ~  0.3  used  based  on  (Sivaramakrishnan 
etal.  2009). 


HDMR  metamodel 


Global  analysis  (750  K) 


UNIVERSITY  OF  LEEDS 


Scatter  plot 


Metamodel 


Wide  distribution  from  zero  to  75( 
Experimental  value  very  infreque 
Missing  steps  or  too  tight  uncertainties  used? 


T/K 


Global  sensitivities:  no  H 

dominant  reaction  UN,VERSITY  of  leeds 


JSR  ( Tres  6  s,  P 1  atm,  stoich,  4%  n- butane,  750  K) 


H202(+M)=0H+0H(+M) 
R26C4H9+02=>C4H8+H02 

ho2+ho2=h2o2+o2  , 
c4h8+oh=>c3h7+hcho 

R41C4H904=0H+C4H803_, 

c4h9=>c2h5+c2h4J 

C4H10+HO2=>H2O2+R20C4H9 

c4h7+ho2=>c4h7ooh 

C4H8+H=>R20C4H9 
R20C4H9+O2=>C4H8+HO2 
R28C  H  0  U=>HO  +C  H  "I 

h2o2+oh=h2o+ho2{ 
c4h10+ho2=>h2o2+c4h9J 

C4H1o+OH=>H2O+R2OC4H0 
HCH0+0H=CH0+H20  j 

hcho+ho2=cho+h2o2J 
ch3o2+h2o2=ch3ooh+ho2  , 
ho2+ch3cho=ch3co+h2o2  - 
ho2+ho2=h2o2+o2  * 

R33C4H902P+02=R41C4H904UP 


Normalised  first-order  contributions  to 
the  overall  variance  of  butane  mole 
fraction. 


t — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — r 

0.00  0.05  0.10  0.15  0.20 


T - T 


Improving  model 
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•  The  ntc  region  appears  to  be  highly  uncertain. 

•  Many  reactions  contribute  to  uncertainty. 

•  Are  there  missing  steps  or  are  assumptions  about  input 
uncertainties  optimistic? 

•  Should  a  hierarchical  Bayesian  approach  be  adopted? 

r-  Start  first  with  H2/02  chemistry. 

>  Then  tackle  HCHO  reactions. 

>  Then  C2  chemistry. 

>  Then  try  to  constrain  reactions  from  C4  scheme. 

•  What  other  experiments  could  help  to  constrain  these 
important  parameters  and  how  do  we  find  out? 


Rate  Coef.  (cm  mol 


H02+H02  (Zhou  et  al.  201 2) 
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T - *l-  -  I  1  I  - T  T  1 

L  2  3  4  5 

1 000/T  (K.) 


TST  calculations. 

Transition  state  ridge 
upward  by  1  kcal/mol  to 
better  fit  experiments  of 
Patrick,  Lightfoot  and 
Kappel. 

Original  fmax  (Baulch  1994) 
was  0.4. 

Probably  could  be  reduced 
using  optimisation 
approaches. 


Aldheydes  +  H02 
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ho2+ch3cho=ch3co+h2o2 

•  Baulch  (2005)  suggests  fmax  of  0.7. 

•  Based  on  a  single  set  of  experimental  measurements  from 
1976  (Colket  et  al.). 

ho2+hcho=cho+h2o2 

•  Baulch  (2005)  suggests  fmax  of  0.4. 

•  Recent  TST  calcs,  from  Li  et  al.  2005 

•  Combined  analysis? 

What  new  sets  of  data  could  help  to 
constrain  these  parameters  better? 


toTr  Vk1 


+  f^x.,) 


2nd  oxygen  addition  reactions  to  form 
02Q00H 
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0.02  k 
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0.01 


/4-factor  for  reaction 
R33C4H902P+02=R41  C4H904UP 


Very  scattered  response 
but  data  hints  that  rate  is 
might  be  lower  than  the 
nominal  value  used  which 
was  based  on  analogy  to 
rate  coefficients  for  propyl 
+  02  from  De  Sain  et  al. 
(2003) 

Clearly  a  key  reaction 
class  driving  fuel  reactivity. 

QOOH+02  studied  for 
propane  and  pentane 
theoretically  by  Bozzelli 
and  Coworkers  but  not  for 
butane. 


2nd  oxygen  addition  reactions  ki 

tO  form  02Q00H  university  of  LEEDS 


•  Rate  used  was  estimated  on  the  basis  of  analogy  to  rate 
coefficients  for  propyl  +  02  from  De  Sain  et  al.  (2003) 

•  Recent  theoretical  calculations  exist  from  Goldsmith  et  al. 
(201 2)  for  the  propyl  +  02  system  suggest  that  for  the  main 
channel  the  reaction  forming  02Q00H  is  slower  than  the 
equivalent  rate  for  propyl  +  02  at  750  K. 

•  An  equivalent  temperature  dependant  study  for  butyl  would 
be  very  useful. 

•  Still  large  uncertainties  for  these  types  of  reactions. 

•  If  basic  C1/H  scheme  was  better  constrained  then  then 
indirect  experiments  might  help  to  constrain  rates  for 
reactions  which  are  difficult  to  study  experimentally. 


Forward  steps 
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•  Could  be  useful  to  collate  information  about  conditions  over 
which  certain  reactions  have  high  sensitivity  coefficients. 

i.e.  parameter  focus  rather  than  experimental. 

•  Would  involve  simulating  a  lot  of  data  sets  over  different 
temperatures  and  pressures  for  a  lot  of  different  fuels  in  a 
consistent  way. 

•  Optimisation  approaches  which  combine  data  sets  are 
clearly  a  way  forward  in  terms  of  reducing  input  uncertainty 
distributions. 

•  Possible  for  H2/02  chemistry  ( last  talk). 

•  How  to  do  it  for  complex  fuels  raises  all  sorts  of  questions. 


USC  University  of 
Southern  California 


1st  Ijrtemational  High-Pressure  Flame  Chemistry  Workshop 


Towards  a  predictive  combustion 
chemistry  model  -  Uncertainty 
propagation  and  minimization 

Hai  Wang  and  David  A.  Sheen 

University  of  Southern  California,  Los 
Angeles,  California,  USA 
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Products 


H20,  C02 
H2,  CO 
C2H2,  soot 


Reaction  Model  Development 

The  Current  Approach 


proliferation  of  models 


log(/c  /  crn^  molecule  1  s1) 


Kinetic  Rate  Parameter  Uncertainties 


H  +  02e0H  +  0  (Rl) 


2000 


1000 


500 


•  Uncertainty  factor  ~1.25 


•  Logarithmic  sensitivity  coefficient 
=  0.24  (ethylene-air,  f-  1,  p  =  1  atm) 


•±5%  (±4  cm/s)  uncertainty  in  predicted 
flame  speed  due  to  Rl  alone 

•  Key  question:  How  do  we  propagate 
uncertainties  in  rate  constants  in 
combustion  simulations? 


05 


103  7  1 1K' 


20 
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Mass  burning  rate  (g  cm 


Uncertainty  Uncertainty  Uncertainty  Uncertainty  Uncertainty  UnCGITtdinty 


Burke,  et  al.  (2010) 
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Uncertainty  Uncertainty  Uncertainty  Uncertainty  Uncertainty  UnCGrtSiflty 


Physical  data 


expt  +  calc 

fundamental 
combustion  expt. 


MUM-PCE 


Method  of  Uncertainty  Minimization  -  Polynomial  Chaos  Expansions 

•  Mathematical  foundation  and  numerical  methods:  Sheen  &  Wang  "Kinetic 
uncertainty  quantification  and  minimization  using  polynomial  chaos  expansions/' 
Combustion  and  Flame,  D0l:10.1016/j.combustflame.2011.05.010. 


•  Model  prediction  presented  as  a  (2-s)  band  of  uncertainty  resulting  from  kinetic 
parameter  uncertainties. 

•  Model  uncertainty  may  be  constrained  by  experimental  data  (ignition  delay, 
species-time  history,  flame  speeds  etc) 


MUM-PCE:  Methods 


•  Stochastic  Spectral  Expansion:  express  kinetic 
parameter  x,  as  a  polynomial  expansion  of  basis 
random  variables 


-  xf] + + 


m 


m  m 


j= i 


&=1  y'=& 


Following  N.  Wiener  (1938),  D.B.  Xiu,  etal.  (2002) 

•  Solution  Mapping:  use  polynomial  response 
surface  to  express  the  relation  between  a 
combustion  response  h  and  x 


N 


N  N 


nr  (x) = n,M + ^ar,x,  +Y^b'jjx,xj 


i— 1 


i= 1  7-* 
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Forward  Uncertainty  Propagation 


> 


n  n  n 

Vr  (x)  =  7Jr,0  +  X  +  X  E  Vi*/  Response  surface  from  solution  mapping 

i= 1  /=!  />/ 


Spectral  representation  of  uncertainty  in  x’s 
(mean  =  0,  s  =  0.5,  each  indep’t  of  others) 


M  M  M 

vr  (x>^) = n,  (x(0))+fx,4 

/=1  /=1  j-i 
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Solution  Mapping  Method 

Fit  a  response  surface  to  the  model 


Central-Composite 
Factorial  Design 


Sensitivity  Analysis 
Based  Design 


10 

G.E.P.  Box,  et  al.  (1978),  Frenklach  et  al.  (1992),  S.G.  Davis  et  al.  (2004) 


Mass  burning  rate  (g  cm 


MUM-PCE  -  Application  in  H2/02  Combustion 


2<7  band 


°0  5  10  15  20  25 


Pressure,  p  (atm) 


•  High-pressure  data  sensitize 
kinetics  of  hydrogen  oxidation. 


•  A  large  number  of  models  outside 
experimental  uncertainty  at  high 
pressures. 

Burke,  et  al.  (2010) 


•  2<j uncertainty  band  calculated  by 
MUM-PCE,  based  on  rate  parameter 
uncertainties. 

•  Models  are  statistical  samples  of 
parameter  uncertainties. 


Sheen  &  Wang  (2011) 


Uncertainty  Uncertainty  Uncertainty  Uncertainty  Uncertainty  UnCGrtSiflty 


Physical  data 


expt  +  calc 

fundamental 
combustion  expt. 


MUM 


rate 

parameter 

uncertainties 


fundamental 
reaction  kinetics 
expt  +  calc 


Physical 

problem 


data 

uncertainty 


fundamental 
combustion  expt. 


Method  of  Uncertainty  Minimization 


x  =  x0  + 

Chemical  model 
+  associated  uncertainty 


N  N  N 

r  (X)  =  ?l  j)  +Xa'-'A'  +SS/,'-'/A>V/ 
i=\  i= 1  j>i 

Physics  model 


(x,5)  =  J]r  (x(0)  )  + 

i=\  i=l  j=i 


z  = 


=1(afs)' 


-(bxo 


*  T -m  *T-m  _  T  *  T  T 

x0  b  +  ax0  b  +  b  x0a  +aa 


)  +  4I 


-1 


Predictions 

+  associated  uncertainty 


a  =  Z1/2 


MUM-PCE  -  Application  in  H2/02  Combustion 

•  Model  uncertainty  constraining 

•  JetSurF  2.0  H2/CO  submodel 

-  14  species,  41  reactions 


Dataset  1:  From  Davis,  et  al.  (2005): 


No. 

P0,  Ps  (atm) 

To,  Ts  (K) 

f 

Laminar  Flame  Speeds 

12 

1-15 

298 

1. 0-3.0 

Ignition  Delay  Times 

13 

0.5-33 

1000-2600 

1.0-6. 1 

Flow  Reactor  Profiles 

9 

1.0-16 

915-1040 

0. 3-1.0 

Laminar  Flame  Profiles 

2 

0.047 

400 

1.9 

Dataset  2: 


From  Burke,  et  al.  (2010): 

No. 

P0,  P5  (atm) 

To,  Ts  (K) 

f 

Laminar  Flame  Speeds 

18 

15-25 

298 

0. 8^-2. 5 

Mass  burning 
rate  (g  cm'2  s'1 


MUM-PCE  -  Application  in  H2/02  Combustion 


Considering  no  experiments 


Dataset  1  (knowledge  prior  to  2010) 


Dataset  1+2  (latest  knowledge) 
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Sheen  &  Wang  (2011) 


Mass  burning  rate  (g  cm-2  s’1) 


MUM-PCE  -  Application  in  H2/02  Combustion 


H2/02/He  mixtures  at 
equivalence  ratio  1 


H2/02/Ar  mixtures  at 
equivalence  ratio  2.5 
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Sheen  &  Wang  (2011) 


MUM-PCE  -  Application  in  H2/02  Combustion 


Dataset  1 
Knowledge  prior  to  2010 


+  Burke,  etal.  (2010) 
Current  knowledge 


H+02<h>0+0H 

h+o2+h2o^ho2+h2o 

h+o2+ar^ho2+ar 

h+o2+he^ho2+he 

H2+02<h>H02+H 

H02+H<h>20H 

oh+ho2^h2o+o2 

co+oh<->co2+h 


* 


Weak  constraint  by  experiments 
Strong  constraint  by  experiments 


H02+H^20H 
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.^\  JetSurF  -  A  Jet  Surrogate  Fuel  Model 


^  Q  ^  Page  -  Safety  - 
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JetSurF  -  A  Jet  Surrogate  Fuel  Model 

JetSurF  is  a  detailed  chemical  reaction  model  for  the  combustion  of  jet-fuel  surrogate.  The 
model  is  being  developed  through  a  multi-university  research  collaboration  and  is  funded  by 
the  Air  Force  Office  of  Scientific  Research.  Project  participants  include 


F.  N.  Egolfopoulos,  HaiWang 

R.  K.  Hanson,  D.  F.  Davidson,  C.  T.  Bowman,  H. 

Pitsch 

C.  K.  Law 

N.  P.  Cernansky,  D.  L.  Miller 
W. Tsang 
R.  P.  Lindstedt 
A.  Violi 


University  of  Southern  California 
Stanford  University 

Princeton  University 
Drexel  University 

National  Institute  of  Standards  and  Technology 
Imperial  College,  London 
University  of  Michigan 


New  Release:  JetSurF  Version  2.0  -  A  working  model  for  the  combustion  of  n-alkane  up  to  n- 

dodecane,  cyclohexane,  and  mono-alkylated  cyclohexane  up  to  n-butyl- 
cyclohexane 

(Release  Date:  September  19,  2010) 


Old  Releases:  JetSurF  Version  1.1  -  A  interim  model  for  the  combustion  of 

n-butyl-,  n-propyl-,  ethyl-,  and  methyl-cyclohexane  and  cyclohexane 
(Release  Date:  September  15,  2009) 


JetSurF  Validation  -  Species  Concentrations 
behind  reflected  shock  waves 


B.  Sirjean,  E.  Dames,  D.  A.  Sheen,  X.-Q.  You,  C.  Sung,  A.  T.  Holley,  F.  N.  Egolfopoulos,  H.  Wang,  S.  S.  Vasu,  D.  F. 
Davidson,  R.  K.  Hanson,  H.  Pitsch,  C.  T.  Bowman,  A.  Kelley,  C.  K.  Law,  W.  Tsang,  N.  P.  Cernansky,  D.  L.  Miller,  A.  Violi,  R. 
P.  Lindstedt,  A  high-temperature  chemical  kinetic  model  of  n-alkane  oxidation,  JetSurF  version  1.0,  September  15, 
2009  (http://melchior.usc.edu/JetSurF/Versionl_0/lndex.html). 


Plot  stolen  from  Ron  Hanson.  Solid  line:  experiments;  dashed  line:  JetSurF 
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Mole  Fraction  Mole  Fraction  Mole  Fraction  Mole  Fraction 


Prediction  Uncertainties  in  As-Compiled  Model 


10-S  10-4  10-3 

Time  (s) 


Good  nominal 
prediction  with 
significant 
uncertainty! 
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Chemistry  Model  &  Experimental  Targets 

•  Modified  JetSurF  1.0 

-  196  species,  1478  reactions 


No. 

P0,  P5  (atm) 

Tv  T5  (K) 

f 

Laminar  Flame  Speeds 

4 

1 

353 

0.8-1. 4 

Ignition  Delay  Times 

11 

1-4 

1000-2600 

0.5-2 

No. 

P5  (atm) 

7-5  (K) 

f 

OH,  H20,  C02,  C2H4,  CH3 
Species  Profiles 

11 

1.6-2. 4 

1365-1545  K 

1 
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Mole  Fraction  Mole  Fraction  Mole  Fraction  Mole  Fraction 


Predictions  of  As-Compiled  and  Uncertainty-Minimized  Models 


Unconstrained 


Constrained 


IQ'5 


10-4 

Time  (s) 


IQ'3 


10-5 


10-4  10-3 

Time  (s) 
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Effect  on  Flame  Speed  Predictions 


0  9  11  1.3 


Considering  no  experiments 


Model  constrained  by  species  profiles 


Model  constrained  by  species  profiles 
+  flame  speeds 


0.7 


Equivalence  Ratio,  <j> 


1.5 
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Effect  on  Flame  Speed  Predictions 


0  5  10  15  20 


1/(2<7obs) 
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What  did  uncertainty  minimization  do? 


Model  constrained  by  species  profiles 


h+o2<->o+oh 

HC0+H^C04H2 
HCO+M++CO+H  +M 
CH-MH2^GH24H 

ch2-hd2<->hco-hdh 

CH2402^C02+2H 
CH  3+H  (+M  )++CH  4(+M ) 
CH340^GH204H 
CHs40H^CH2-4M20 
2GH  3(+M  )++C2H  6  (+M ) 
c2H3-m^c2H2-m2 

C2H 3402^GH  2GH  O  -HD 
C2H4-HD^CH3-HHCO 
C2H4-hDH^G2H3-HH20 
aC^H  g+H  (+M  )++  C^H  g  (+WI ) 
C^H  6  +H  <-4aCgH  g+H  2 

PXG^11^2H5^NC7H16 
pG4H  9 +nC^H  ?++N  C?H  1 6 

NG?H  16-HH<->PXC7H  1  g+H  2 
NC?H  16-HH++SXC7H 1  g+H  2 
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Model  constrained  by  flame  speeds 


CH3,  CH2,  secondary  chain  H  chain  branching 

branching,  fuel  breakup 
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What  did  uncertainty  minimization  do? 
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"Our"  Approach 


New 

Sensitivity 

experiments 

analysis 

Electronic 

structure 

calculations 

Reaction 
rate  theory 


Uncertainty 

minimization 


Response  surface 
development 


Take  home  messages: 

•  Consistency  is  more  critical 
than  predictability! 

•  Need  a  comprehensive 
model  to  truly  reflect  the 
inherent  model  hierarchy 
and  uncertainty. 
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VELAAft 


Measurements  of  Elementary  Rate  Constants, 
Ignition  Delays  and  Species  Histories  in  Shock  Tubes 


R.  K.  Hanson,  D.  F.  Davidson 
Stanford  University 

1st  International  Workshop  on  Flame  Chemistry 
July  28-29,  2012 


•  Shock  Tube/Laser  Approach 

•  Advances  in  Methodology 

•  Elementary  Reaction  Rate  Studies 

•  Multi-Species  Time-Histories 


Kinetics 
Shock 
Tube  1 
(30  atm) 


Kinetics 
Shock 
Tube  2 
(30  atm) 


Stanford  Shock  Tubes 


Aerosol 
Shock 
Tube 
(10  atm) 


High 
Pressure 
Shock 
Tube 
(500  atm) 


New  Constrained  Reaction  Volume  (CRV)  Facility 

2  Configurations 


Sliding 

Valve 


Constrained 
Reaction  Volume 


CRV-1 

Gas-Phase 
Experiments 
to  100  atm 


CRV-2 

Aerosol 
Experiments 
to  100  atm 


Aerosol 

Generation  Tank 
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First  use  of 
tunable 
dye  lasers 
in  shock 
tubes 
(1982) 


Ultra-fast 
lasers  used  to 
extend  UV 
tuning  range 
(2009) 


Stanford  Laser  Diagnostics 


New  lasers 
allow  simple 
access  to  mid-IR 


Ultraviolet  Visible  Infrared 


ch3 

216  nm 

CN 

388  nm 

H20 

2.5  jum 

NO 

225  nm 

CH 

431  nm 

C02 

2.7  jam 

o2 

227  nm 

NCO 

440  nm 

ch4 

3.4  |iim 

ho2 

230  nm 

no2 

472  nm 

ch2o 

3.4  |iim 

OH 

306  nm 

nh2 

597  nm 

CO 

4.6  |iim 

NH 

336  nm 

HCO 

614  nm 

NO 

5.2  |iim 

c2h4 

10.5  |um 

(2007-10) 


Detection  Limit  [ppm] 


Laser  Absorption  Yields  High  Sensitivity 

Representative  Detection  Limits:  Polyatomic  Molecules 

1000 

100 

10 

1 

0.1 

0.01 

1000  1500  2000  2500  3000 

Temperature  [K] 

*  Polyatomic  molecules  @  1500K:  2-200  ppm 


Laser  Absorption  Yields  High  Sensitivity 

Representative  Detection  Limits:  Diatomic  Molecules 


sub-ppm 
sensitivity  for 
OH,  CH,  CN 


Diatomic  molecules  @  1500K: 

-  sub-ppm  detectivitiy  for  UV  absorbers 

-  ppm  detectivity  for  IR  absorbers 
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Advances  in  Shock  Tube 
Methodology 

1.  Improved  uniformity  with  driver  inserts 

2.  Longer  test  times  with  tailored  gas 
mixtures  &  extended  driver 

3.  Reactive  gas  modeling: 
problem  and  solutions 


Improvement  in  Reflected  Shock  Temperature  Uniformity 

Using  Driver  Inserts 


Conventional  shock  tube 
operation  can  provide 
near-ideal  uniform  flows 
for  1-3  ms 

But,  boundary-layers  and 
attenuation  induce  dP/dt 
and  dT/dt  at  longer  times 


0.8%  Propane/  8%  OJ  Ar 

J _ i _ I _ i _ I _ i _ I 

0  5  10  15 

Time  [ms] 
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Improvement  in  Reflected  Shock  Temperature  Uniformity 

Using  Driver  Inserts 


These  effects  reduce 
ignition  delay  times 
relative  to  Constant  P  case 

Solution:  Driver  Inserts 


Driver  inserts  modify  flow  to  achieve 
uniform  T  and  P  at  long  test  times 
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Improvement  in  Reflected  Shock  Temperature  Uniformity 

Using  Driver  Inserts 


Result:  dP/dt  =  0  prior  to 
ignition 

Proper  xign  for  comparison 
with  Constant  P  simulations 


1 

dPIdt- 

l. 

ign 

=  0%/ms  1 

Li 

^initial  =  1039  K 

0.8%  Propane/  8%  OJ  Ar 

1 . l . l_ i_ l 

0  5  10  15 


Time  [ms] 
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Longer  Test  Times  Achievable  with 
Tailored  Gas  Mixtures  &  Extended  Driver  Sections 


10 


*  Conventional  shock  tube  operation:  ~  1-3  ms  test  time 

*  No  overlap  with  RCM  operation~  10-150  ms  test  time 


0) 

IE-3  I - 1 - 1 - 1 - J- 

0.4  0.8  1.2  1.6 

1000/T  [1/K] 


ii 


Longer  Test  Times  Achievable  with 
Tailored  Gas  Mixtures  &  Extended  Driver  Sections 

•  Longer  driver  length  and  tailored  gas  mixtures 
can  provide  longer  test  times  (>  40  ms) 


2x  Driver  Extension 


1000 


2500K 


1000K 


625K 


“  100  - 


10  - 


0.01  - 


IE-3 


0.8  1.2 
1000/T  [1/K] 


Shock  tubes  now  can  overlap  with  RCMs 
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First  Use  of  Long  Test-Time  Facility: 

Low  Pressure  n-Heptane  Ignition  in  the  NTC  regime 


6 

Test  time  =  45  ms  at  725K 

5 

Enables  first  low-pressure  § 

(~3  atm)  n-heptane  ignition  § 

data  in  NTC  regime  3 

> 

_C0 

0)  9 

Clear  evidence  of  2-stage 
ignition:  xx  =  8  ms,  x2  =  20  ms 

Next  step:  add  driver  inserts  to  o 
remove  dP/dt  &  c/T/c/t 


How  do  the  3  atm  data 
compare  to  high  P  data? 


Ignition  Delay  Times  [ms] 


Low-Pressure,  Low-Temperature  Studies: 
New  Results  for  n-Heptane  in  the  NTC  Regime 


Previous  12-55  atm  data 
New  3  atm  data 

How  do  measurements  compare  with  current  models? 


Low-Pressure,  Low-Temperature  Studies: 
NTC  Heptane,  Comparison  with  Model 


•  LLNL  C7  (2000)  model  performs  reasonably  well  at  high  P 

•  Further  tests  needed  at  low  P 

•  Reveals  value  of  long  test  time  experiments 
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Reactive  Gasdynamics  Modeling:  A  Problem 


*  Most  current  reflected  shock  modeling  assumes 
Constant-Volume  or  Constant-Pressure 


But: 

•  Exothermic  energy  release  during  oxidation  or 
endothermic  cooling  during  pyrolysis  changes  T  &  P 
behind  reflected  shocks 

not  a  Constant-V  or  Constant-P  process! 

•  Example:  Heptane  Ignition 
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Effect  of  Energy  Release  on  P  Profiles: 
n-Heptane  Oxidation 


Measured 
Sidewall  P 

Constant  P 
Model 


-100  0  100  200  300  400  500  600 

Time  [us] 


How  does  this  compare  with  models? 
Not  a  constant  P  process! 
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Effect  of  Energy  Release  on  P  Profiles: 
n-Heptane  Oxidation 

Constant  V 
Model 

Measured 
Sidewall  P 

Constant  P 
Model 


-100  0  100  200  300  400  500  600 

Time  [us] 


How  does  this  compare  with  model? 

Not  a  constant  P  process! 

Not  a  Constant-V  process,  even  for  0.4%  fuel! 
So  how  can  entire  process  be  modeled? 


3  Proposed  Solutions  to  Enable  Modeling  through 

Entire  Combustion  Event 

1.  Minimize  fuel  loading  to  reduce  exothermically-  or 
endothermically-driven  T  and  P  changes 

-  enabled  by  high-sensitivity  laser  diagnostics 

2.  Modified  gasdynamics  modeling  to  account  for 
P  and  T  change  during  combustion 

-  work  in  progress  (but  computationally  intensive:  1-D,  3-D) 

3.  Use  new  constrained  reaction  volume  concept  to 
minimize  pressure  perturbations 

-  enables  constant  P  (or  specified  P)  modeling 

Examples:  1)  Use  of  dilute  reactive  mixtures 

2)  Use  of  constrained  reaction  volume 
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Pressure  [atm] 


Example  1:  Benefit  of  Dilute  Mixtures 
3-Pentanone  Oxidation 


Low  Fuel  Loading  Experiment 


400  ppm  3-Pentanone/02/Ar 
1486  K,  1.52  atm,  <[>=1 

_l _ i _ I _ i _ I _ i _ 

0  500  1000  1500 

Time  [|as] 


OH  Mole  Fraction 


Time  [|is] 


*  Pressure  nearly  constant  throughout  experiment 

*  Good  agreement  between  Constant  H,P  model  and  expt. 

*  Model  successfully  includes  temperature  change 
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Example  2:  Constrained  Reaction  Volume  Approach 

Hydrogen  Ignition  at  950  K 


Conventional  Shock  Tube 

Pre-Shock 


Reflected 
Shock  Wave 


Constrained  Reaction  Volume 

Pre-Shock 


Helium 

Test  Mixture 

Helium 

Non-Reactive  Mix 

TM 

Post-Shock 

Post-Shock 

—ym 

*1 

i 

Large  Region  of 
Energy  Release 


Small  Region  of 
Energy  Release 


*  Large  reaction  volume  gives  large  energy  release  AP  &  AT 

*  CRV  gives  reduced  energy  release  near-constant  P 
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Example  2:  Constrained  Reaction  Volume  Approach 

Hydrogen  Ignition  at  950  K 


Conventional  Shock  Tube  Constrained  Reaction  Volume 


*  Conventional  ST  exhibits  large  pressure  change! 

*  CRV  pressure  nearly  constant  throughout  experiment! 

*  Allows  kinetics  modeling  through  ignition  and  combustion! 
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Elementary  Reaction  Rate 
Determinations 


Goal 

Near-direct  determination  of  rate 
constants  for  specific  reactions 

Examples 
OH  +  ketones: 

acetone,  2-butanone,  2-  &  3-pentanone 
OH  +  alkanes 
OH  +  butanol  isomers 
OH  +  methyl  esters 


Experimental  Strategy 


TBHP  used  as  a  prompt  OH  precursor 

*  Useful  T  range  (850  to  1350  K) 

*  Pioneered  by  Bott  and  Cohen  (1984) 

*  Also  used  at  Argonne 

Fuel  in  excess,  pseudo-first  order  experiment 


CH3 

I  / 

H.C-C-0 

I 

CH, 


OH 


fast  upon 
shock  heating 


o 

11 

h3c  ch3 


ch 


tert-buty  I  hydro  peroxide 


Acetone 


o 


Representative  3-Pentanone+OH  Data: 
1188  K  and  1.94  atm 


•  High-quality  data  allows  high-precision  comparison  with  model 
9  Sensitivity  analysis  confirms  pseudo-first  order  behavior 
9  Near-direct  determination  of  reaction  rate  constant! 
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Results  for  3-Pentanone  +  OH  Products 


No  other  high  temperature  data  available 

NUI  model  (Serinyel  et  al. )  in  excellent  agreement 


Summary:  OH+Ketones-^ Products 


1429K  1250K  1 1 1 1K  1000K  909K  833K 


1000/T  [1/K] 

How  do  data  compare  with  Structural  Activity  Relationship  (SAR)  model? 
Data  agree  within  25%  with  SAR-estimated  rate  constants 
Similar  measurements  performed  with  methyl  esters 


Methyl  Ester  +  OH 


Summary:  OH+Methyl  Esters  Products 


1429K  1250K  1 1 1 1K  1000K  909K  833K 


Data  agree  within  25%  with  SAR-estimated  rate  constants 
Current  work:  OH  +  aldehydes,  alcohols 


Multi-Species  Time-Histories 


Motivation 

Multi-species  provide  greater  constraint 
on  mechanism  refinement/evaluation 

Recent  Work 

Ketones:  Acetone,  Butanone,  3-Pentanone 
Alcohols:  1-,  2-,  tert-,&  iso-Butanol 
Alkanes:  n-Hexadecane 


Esters:  Methyl  Formate,  MA,  MP,  MB,  EP 


3-Pentanone  Mole  Fraction  [%] 


Multi-Species  Approach:  3-Pentanone  Pyrolysis 


Excellent  SNR,  high  sensitivity  data 
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3-Pentanone  Mole  Fraction  [%] 


Comparison  with  Serinyel  et  al.  (2010)  Galway  NUI 


Two  Differences:  1)  3-P  decomposition  rate;  2)  CO/C2H4  yields 
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3-P  Data  Enables  Revision  of  Decomposition  Rate 

3-Pentanone  Sensitivity  Arrhenius  Plot:  3-Pent  Products 


1429K  1250K  1 1 1 1K 


3-P  data  show  strong  sensitivities 
to  k,  +  k 


Revised  ktota ,  3.5x  Serinyel  et  al.  rate 
CO  yields  still  not  correct! 
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Mole  Fraction  [%] 


CO  Yield  Resolved  through  Use  of  O-Atom  Balance 


3-P  and  CO  data  yield  total  O-atoms 


O-atom  concentration  at  1.5  ms 

Laser  Absorption 

3-Pentanone:  23% 

CO:  69% 

Sum:  92% 

Simulation  (with  new  k^k2) 

3-Pentanone:  23% 

CO:  43% 

CH3CHCO:  27% 

Sum:  93% 


Model  underpredicts  CO  and  overpredicts  methyl  ketene 

Why?  Methyl  ketene  decomposition  pathway  missing  in  mechanism 

Introduce  CH3CHCO  ->  C2H4  +  CO  (assume  k  the  same  as  ketene  decomp.) 
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CO  Mole  Fraction  [ppm] 


Revised  Model  Improves  3-P  Simulations 


Final  Modifications  to  Serinyel  et  al. 
3-Pentanone  Mechanism 

*  Revised  decomposition  rate: 
3-pentanone  products 

*  Additional  reaction: 
methylketene  C2H4  +  CO 


CO,  and  Low-T  C2H4  time-histories 
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Ongoing  Work 


Diagnostics  development  &  spectroscopy: 

-  aldehydes  (CH20,  CH3CHO) 

-  alkyl  radicals  (C2H5) 

-  methyl  esters  (methyl  formate) 

-  Alkenes  (C3H6,  C4H8) 

Continued  improvement  of  shock  tube  methods: 

-  constrained  reaction  volume 

Direct  measurement  of  elementary  reactions: 

-  OH  +  oxygenates  (aldehydes,  ethers,  alcohols) 

-  decomposition  of  oxygenates 

-  CH3,  H02  reactions  with  HC 
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crf  The  role  of  advanced  diagnostics  in 
combustion  chemistry  research 
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Introduction:  The  problem 


i 


m 


>Clean  combustion:  high  efficiency,  low  emissions 
Alternative  transportation  fuels:  new  chemistry 
>Novel  combustion  regimes:  influences  of  p,  T,  phi,  mixture 


Needed:  Systematic 
knowledge  on 
combustion  chemistry 


o 

03 


a> 

o 

c 

ca 

> 


v  3 


cr  2 
LU 


1  • 


O' 
100 


140  180  220  260  300 

Temperature 
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Probing  flame  chemistry:  Wanted 


i 


m 


AX 

sfesfe 


n-C4H10 


iso-C4H10 


>ln  situ  information  on 
unknown  species  mix: 
identification. 

>  Large  dynamic  range  of 
mole  fractions,  labile 
species. 

>Quantitative  set  of 
concentration-reaction 
time  profiles. 

>  Example:  butane 
combustion,  ~  40  species. 


66% 


49% 


OfBwald  et  al,  Z  Phys  Chem  225,  201 1 ,1029 


100% 


33% 
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C2,  cyclopentene 
■  C/O  =  0.5 
a  C/O  =  0.6 

*  C/O  =  0.7 

•  C/O  =  0.77 


C2,  propene 

• 

C/O  =  0.4 

■ 

C/O  =  0.5 

▲ 

C/O  =  0.6 

* 

C/O  =  0.7 

• 

C/O  =  0.77 

CH, 

propene 

• 

C/0  = 

0.4 

■ 

C/0  = 

0.5 

▲ 

C/0  = 

0.6 

* 

C/0  = 

0.7 

> Laser  spectroscopy:  Raman,  LIF,  CRDS,  IR  absorption,  etc. 
>Quantitative,  non-intrusive,  small  molecules:  e.g.  CH,  C2 
>Not  suitable  for  complete  flame  analysis 
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0  2  4  6  8  10  12  14 

h  /  mm 


M.  Kohler  et  al,  J  Phys  Chem  A1 14,  2010,  4719 
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Probing  flame  chemistry:  laser  diagnostics 


ll  f /, 


1st  International  Workshop  on  Flame  Chemistry,  Warsaw,  Poland,  July  28-29,  2012 


4 


qdd  /  X  qdd  /  X 


Universitat  Bielefeld 

Department  of  Chemistry 

.Ml 

m 

* 

r 

Small  species  and  chemical  mechanisms 

\  _ 

1st  International  Workshop  on  Flame  Chemistry,  Warsaw,  Poland,  July  28-29,  2012 
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Universitat  Bielefeld 


Probing  flame  chemistry:  mass  spectrometry 


> Molecular-beam  mass  spectrometry  (MBMS): 

>  Complete  species  set  (?!);  structure-sensitive,  isomers 
>Quantitative  -  but:  fragmentation,  overlaps,  sampling . 


Distance  from  Burner 


I  urbo 
pump 


W1UP  .. 

TQF  tube 


Ll_  C_ 


Rerlectron 


Pressure 

control 


^  Turtib 

a  pump 

E 


Flat 

bunsr 
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Ionization  chamber 


Relative  Photoionization  Efficiency 


Department  of  Chemistry 

Species-  < 

and  structure-selective  analysis 

wm  a 

A  *&&& 

>Tunable  VUV  single-photon  PI-MBMS: 

> Distinction  of  intermediates  in  C/H,  C/H/O,  and  C/H/O/N  systems 
>Mass  &  photoionization  efficiency  (PIE)  spectra  as  discriminators 


9.0  9.5  10.0  10.5  11.0 

Photon  Energy  (eV) 

C.A.  Taatjes  et  al,  PCCP  10,  2008, 1 


8.0  8.5  9.0  9.5  10.0 


Photon  Energy  /  eV 

B.  Yang  et  al,  Combust  Flame  148,  2007,198 


8  9  10  11  12 


Photon  Energy  (eV) 

m/z=41 ,  C3H5,  C2H2NH,  CH3NC  (CH3CN) 
A.  Lucassen  et  al,  PROCI  34,  2012 


1st  International  Workshop  on  Flame  Chemistry,  Warsaw,  Poland,  July  28-29,  2012 
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M 

Combining  high  mass  +  energy  resolution 

wm  a 

A  *&&& 

>Fuel  complexity: 
hydrocarbons  — ► 
oxygenates  — ►  fuel-N. 

>More  elements  need 
better  mass  separation: 
— ^EI-MBMS. 


20  40  60  80  100  8.5  9.0  95  10.0  10.5  11.0  11.5  12.0  12.5  13.0  13.5  14.0 

m/z  photon  energy  /  eV 


>Fuel  structure  and 
intermediate  mix  need 
isomer  separation: 

— >PI-MBMS. 

>  Example:  morpholine 
flame. 


A.  Lucassen  et  al,  PROCI  32,  2009,1269 


1st  International  Workshop  on  Flame  Chemistry,  Warsaw,  Poland,  July  28-29,  2012 
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mole  fraction 


Universitat  Bielefeld 

Department  of  Chemistry 

mm 

▲  .  * 

'/mi 

'  -4 

Combining  EI-MBMS+PI-MBMS+LIF+CRDS 


>MBMS  results  for  HCN,  NH3,  NO,  HNCO,  CH3, ... 
>LIF  temperature  measurement 
>CRDS  results  for  NH2,  CN,  CH,  OH,  ... 


h  /  mm 


h  /  mm 

P.  Nau  et  al,  Exp  Fluids  32,  2009,1269 


1st  International  Workshop  on  Flame  Chemistry,  Warsaw,  Poland,  July  28-29,  2012 
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Experiment  vs.  model:  (un)disturbed  T  profile 


_ 

Bfrdi  *  * 


>  T profile  for  optical  or 
MBMS  sampling  must  be 
considered  in  comparison 
with  models  -  no  shifts! 

>  Example:  butanol  flames. 


HAB  (mm) 


M.  Sarathy  et  al,  Combust  Flame  159,  2012,  2028 


1st  International  Workshop  on  Flame  Chemistry,  Warsaw,  Poland,  July  28-29,  2012 
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Temperature  (K) 


Universitat  Bielefeld 


Department  of  Chemistry 


New  combustion  details:  study  of  isomeric  fuels 


mm 

5£'4fiS 


^Quantitative  but 
relative  species 
measurement  is 
more  forgiving. 

>Chemical  trends 
are  more  obvious. 

>Model  can  probe 

mechanistic 

differences. 

>  Example:  ethanol 
and  DME  addition 
to  propene. 


A.  Frassoldati  et  al,  Combust  Flame  158,  2011,  1264 


1st  International  Workshop  on  Flame  Chemistry,  Warsaw,  Poland,  July  28-29,  2012 
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Fuel-N  conversion:  biomass  components 
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HC-NH-CK 
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>  Isomeric  fuel 
decomposition 
pathways  different. 

>Ethylamine  break¬ 
down  suggests  early 
NH3  formation. 

Substantial  HCN  mole 
fractions  expected  for 
both  fuels. 

>  Example:  ethylamine 
and  dimethylamine 
combustion. 


A.  Lucassen  et  al,  Combust  Flame  159,  2012,  2254 
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Ethylamine  vs.  dimethylamine  combustion 
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A.  Lucassen  et  al,  Combust  Flame  159,  2012,  2254 


>NH3  mole 
fractions  much 
higher  for  EA. 

^  Model  under¬ 
predicts  NH3  in 
both  flame  sets. 

>High  HCN  mole 
fractions  of  up 
to  1 2%. 

>Model  under¬ 
predicts  HCN. 
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Low-temperature  combustion  chemistry 

\ 

wm  a 

(a) 

To  pump 

Skimmer 

Quartz  nozzle 
Heating  wire 


Spiral  driving 
localizer 


To  TOF 


Step  motor 


(C) 


Premixing  zone 


Reaction  zone 


CO2 


> Sampling  from  highly-diluted,  preheated,  partially-premixed 
DME  flame  at  1  bar  with  T  =  1400  K. 

>Partially-premixing  species  profiles  are  position-dependent. 

K.  Zhang  et  al,  PROCI  34,  201 2  I 
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Mole  Fraction  (xIO'3) 


Department  of  Chemistry 


K.  Zhang  et  al,  PROCI  34,  2012 
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Low-temperature  combustion  chemistry 


>Highly  diluted 
preheated  DME 
combustion. 

>Three  stoichio¬ 
metries:  0.85, 
1.0,  and  1.2. 
Intermediate-T 
behavior;  e.g. 
CH3  vs.  CH20. 


1 E06  Mon 


DME 


(a)  <D  =  0.85 


(d)  0  =  0.85,  CH20 


(g)  O  =  0.85,  CH3 


(b)  O  =  1.00 


DME 


(e)  O  =  1.00,  ChLO 


K  (h)  O  =  1.00,  ch3 


(c)  O  =  1.20 


(f)  0  =  1.20,  CH20. 


(i)  O  =  1.20,  CH 


-a—  Position  1 
3— Position  2 
Position  3 
Position  4 
Position  5 
Position  6 
Position  7 
Position  8 
Position  9 

0.0  05  To  15  2^0 


Height  Above  Burner  (mm) 
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Low-temperature  combustion  chemistry 


% 


180  mm 

I  I 


heating  zone:  740  mm 


non-heated  zone:  540  mm  40  mm 

I  I 


UXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX1IXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX1 


fuel, 


02,  Ar 


F=> 


inner  diameter 
8  mm 


•T  ◄- 

gas  ^ 


350  mm 


reactor  length:  1500  mm 


H  exhaust  gas 


CO2:  O  •  C  O2:  A  A/\  CO:  □  ■  B 


F.  Herrmann  et  al,  PROCI  34,  2012 


>Highly  diluted 
DME  oxidation 
in  flow  reactor. 

>Three  stoichio¬ 
metries:  0.8, 1.0, 
and  1.2. 

Expected  low-T 
behavior. 

>Compare  with 
model/  to  EtOH: 
— >WIP! 

1E07  Mon,  W5P070  Fri 
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Low-temperature  combustion  chemistry 


m/z 


>DME  oxidation  in 
low-T  flow  reactor. 

>High  mass  reso¬ 
lution:  fuel  DME 
(46.042)  and  formic 
acid  (46.01 0)  at  Th  = 
583  K  are  separated. 

>  Methyl  formate  is 
detected. 

>More  species,  VUV- 
PI-MBMS:  — >WIP! 

I — N  1 E07  Mon,  W5P072  Fri 
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M.  Schenk  et  al,  Combust  _ 

Flame  201 2,  submitted  = 

o> 
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Sampling  at  p  =  10_1-10'6  mbar 


>ln-situ  isomer  separation  +  EI-MBMS 
>  Example:  butene  flames,  C5H8  m/z  68 


Quartz  nozzle 
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Butene  combustion  as  C4  chemistry  subset 


% 
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>  Different  fuel  break¬ 
down  schemes  for  the 
3  butene  isomers. 

>C3  route  is  of  high 
importance  for  all 
butenes,  almost 
exclusive  for  /-butene. 

> Additional  C4  route 
exists  for  the  2  linear 
butenes. 


W2P081  Tue 


M.  Schenk  et  al,  Combust  Flame  2012,  submitted 
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Highly  complex  chemistry:  Biomass  pyrolysis 

>Tunable  VUV  PI-MBMS  for  analysis  of  pyrolysis 
profiles  (T,  time)  of  second-generation  biofuels. 


>Fast-growing  poplar  wood  as  potential  corn 
replacement;  carbon  conversion  to  biofuels. 


J.  Weng  et  al,  PROCI  34,  2012 
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1 

1 

Poplar  pyrolysis  set-up 

m 

wm  a 

1 

>Qstar  triple  quadrupole 
TOF-MS  with  mass  range  30- 
20000  Da,  resolution  10000. 

>  Reactor  is  heated  to  specific 
temperature,  then  sample 
inside  the  quartz  pole  is 
pushed  into  the  furnace. 

>  Pyrolysis  products  pass 
through  a  repeller  plate  into 
the  photoionization  region, 
photoions  are  analyzed  by 
the  QTOF  mass  spectrometer. 


II.  Photonionization 

J.  Weng  et  al,  PROCI  34,  2012 
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C4H602  (86)  /  C6H60  (94)  H  C5H10O2  (102)  /  C7HgO  (108)  H  C6H602  (110)  H  C5H603  (114)  /  C7H802  (124)  G 

2,3-Butanedione  Phenol  Valeric  acid  o-Cresol  Dihydroxybenzenes  4-Hydroxy-5, 6-dihydro  2-Methoxyphenol 

-2H-pyran-2-one 


OH 


C6H603  (126)/ 
5-Hydroxymethyl- 
furfural 


"CHO 
OH 

C8H802  (136)  H 
2-Hydroxy-5- 
methylbenzaldehyde 


vOCH3 
OH 
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4-Methylguaiacol 


OCH, 
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4-Vinylguaiacol 


OH 
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CHO 

C8H803  (152)  G 
Vanillin 


OCH, 


C9H1202  (152)  G 
4-Ethlguaiacol 


H3CO  y  och3 
OH 

C8H10O3  (154)  S 
2,6-Dimethoxyphenol 


OH 


OCH. 


9  X 


OH 


OCH3 

W  CoH 


OH 

H3CO  och3 


och3 
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OCH, 


OH 


h3co 


och3 


C6Hi0O5  (162)/ 
Levoglucosan 


Ci0H12O2(164)G 

Eugenol 


Ci0H12O2  (164)  G  C9H10O3  (166)  G  C9H1203  (168)  S 
Isoeugenol  Apocynin  4-Methyl-2,6- 

dimethoxyphenol 


CioHi0O3  (178)  G 
Coniferyl  aldehyde 


Ci0H12O3(180)S 

Vinylsyringol 


.OH 


OCH3 


Ci0Hi2O3  (180)  G 
Coniferyl  alcohol 


H3CO  y  OCH3 
OH 

C9H10O4  (182)  S 
Syringaldehyde 


CuHu03  (194)  S 
4-Propenylsyringol 


H3CO  y  OCH3 
OH 


h3co  y  och3 

OH 


C10H12O4  (196)  S  CnH1204  (208)  S 

1-(4-Hydroxy-3,5-dimethoxy  Sinapyl  aldehyde 
-phenyl)ethanone 


J.  Weng  et  al,  PROCI  34,  2012 


Ci1H1404  (210)  S 
Sinapyl  alcohol 
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Time-resolved  pyrolysis  profiles 
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MBMS  in  non-premixed  flames 

wm  a 

Av': 

**. 

9k 

S.A.Skeen  et  al,  PROCI  34,  2012 


Ni  skimmer 
&  quartz 
probe/cone 


Angular 

translator 


VUV 

photons 


Horizontal 

translator 


>C2H2  flame  at  30  Torr 
>Mole  fractions  vs.  fuel  outlet 
>  Agreement  with  model  by  Jim  Miller 
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MBMS  in  non-premixed  flames 
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>  Radicals  can  be  detected. 
>Profiles  appear  shifted  vs.  model. 


5E07  Fr 

S.A.Skeen  et  al,  PROCI  34,  2012 
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Mame-sampiing  aerosol  mass  spectrometry 
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Flame-sampling  aerosol  mass  spectrometry 

£ 

>Gas  phase  species  and  50-100  nm  particles  are  sampled 
from  700  Torr  counterflow  flame  with  quartz  microprobe. 


focused  with  aero¬ 


dynamic  lens  onto 
heated  copper  plate 
and  flash-vaporized. 

>Molecular  con- 


Ion  Optics 


>  Particles  are 


stituents  are  VUV- 
photoionized  and 
detected  by  TOF 
mass  spectrometry. 


1st  International  Workshop  on  Flame  Chemistry,  Warsaw,  Poland,  July  28-29,  2012 


27 


Universitat  Bielefeld 

Department  of  Chemistry 

M 

Particle  chemical  composition:  surprises 

wm  a 
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>  Mass  spectra  map  chemical  composition  of  the  particles  vs. 
distance  from  the  fuel  outlet;  PIE  spectra  permit  identification. 

> Example  m/z=116,  normally  only  identified  as  indene,  but  PIE 
curve  shows  also  phenyl-substituted  allene  and  propynes. 
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Soot  formation  beyond  benzene 
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>Many  species  beyond 
benzene  to  ~1000  Da. 

>Mass  spectra  peak  at 
around  202  Da,  i.e. 
pyrene. 

>Detailed  analysis  in 
progress. 
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Perspectives:  Future  trends  and  needs? 


>New  combustion  chemistry:  more  elements,  higher 
mass  range,  structure-selective  quantitative  analysis. 

>New  measurement  domains:  larger  p  and  T  range, 
pyrolysis  vs.  highly-diluted  oxidative  systems, 
particles,  aerosols,  homogeneous  vs.  heterogeneous 
combustion  chemistry. 


>  Time-resolved  analysis. 


>  Support  by  theory:  ionization  energies,  structures, 
simulation  of  spectra,  kinetics,  thermochemistry. 

>  Combination  of  methods:  optical&MS  techniques,  in 
situ  GC-MBMS,  IR-absorption,  MS-MS,  PEPICO,  .... 


> Caveat:  Let‘s  get  sampling  problem  solved! 
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Elementary  Kinetics  Discussion  Summary 

Robert  Tranter 

•  A  number  of  topics  were  discussed  extensively.  The  main  points  are 
sumarized  on  the  following  slides  and  some  principle  points  are  given 
below.  Interestingly,  the  group  independently  arrived  at  many  conclusions 
similar  to  those  of  the  Reaction  Mechanism  discussion  group. 

•  A  primary  topic  of  discussion  concerned  making  experimental  and 
theoretical  data  available  in  a  format  that  would  enable  future 
reevaluation. 

•  A  secondary  topic  of  some  concern  was  the  long  term  status  of  the  NIST 
kinetics  database.  All  participants  in  the  discussion  acknowledged  the 
importance  of  this  database  to  their  research  and  it  is  clearly  a  significant 
tool  for  the  kinetics  research  community. 

•  Concern  was  expressed,  during  presentation  of  the  discussion  groups 
findings,  that  there  is  only  a  very  limited  and  shrinking  number  of  people 
measuring  elementary  rate  coefficients  at  flame  temperatures. 


•  Sharing  Data 

-  Common  structure 

•  XML  type? 

-  Infrastructure 


PRIME  type? 

More  use  of  supplerru 


-  Attach  XML  files 

—  ASSCII  files  of  data 


reporting  data 


information  in  journals 


NIST  Kinetics  database 

-  What  is  the  current  status? 

-  What  will  be  the  long  term  availability? 

-  Can  the  Combustion  community  provide  support? 

-  Can  it  made  more  open  to  help  keep  it  up  to  date? 

•  Upload  your  data 

•  Feedback  forms 

•  Wikipedia  type  approach? 


Areas  of  need  for  flame  chemistry 

-  Low  temperature  chemistry 

-  H02  reactions 

•  Experiment 

•  Theory 

-  R  +  02 

•  Experimental  measurements 


•  Core  reaction  mechanism 

-  CO,  Cl,  C2.... 

-  Can  a  definitive  set  be  established? 

-  Is  a  single  core  mechanism  desirable? 

-  How  would  it  be  maintained  and  updated? 


Kinetic  mechanisms"  discussion 

Summary 


Matt  Oehlschlaeger 


Challenges 

Expense  of  theoretical  and  experimental 
determinations  of  rate  coefficient 

-  ~30  rate  theoretical  rate  predictions  /  year; 
similar  annual  experimental  yield 

-  Not  all  theoretical  predictions  are  of  equal  value; 
QC  tools  in  the  wrong  hands  can  be  dangerous 

Development  and  optimization  of  rate  rules 
for  reactions  that  can't  be  calculated  or 
measured 

-  Is  number  of  rules  finite? 


Challenges 

Never  will  be  enough  data  to  validate  a  model  across  all 
conditions  found  in  a  turbulent  flame 

-  Need  to  identify  regimes  where  uncertainty  &  sensitivity  exists 

-  Multiple  groups  perform  complimentary  experiments 

-  (Re)evaluate  old  data 

At  high  temperatures 

-  Pyroltic  chemistry  leaves  fragments 

-  Structure  of  fuel  fragments  (olefins)  controls  their  reactivity 

-  Design  thoughtful  experiment  to  probe  these  portions  of  the 
models;  not  always  necessary  to  look  at  parent  fuel 

Low  temperatures:  fuel  specific  chemistry  important 

C0  chemistry 

-  Relatively  small  uncertainties  (~25%)  ->  high  cost  for  small 
improvements 

-  H02  chemistry  a  question,  measurements  needed 


Challenges  -  What  does  industry  care  about? 

Real  fuels 

-  Blending  effects  -  classes  of  components 

-  Prediction  of  octane  rating  for  fuel  blends 

-  Aromatics,  olefins 

-  N-alkanes  (NTC  energy  release) 

-  Alternative  fuels  farther  off  than  public  may  think 

Mechanisms 

-  Simplicity  -  number  of  reactions 

-  Rate  of  heat  release,  transition  from  low-  to  high-T,  temperature  distribution, 
NOx,  soot,  noise  (diesel),  ignition  delay 

-  Resolve  trends 

-  Sensitivity  of  parameters  on  operation/performance 

Tools  for:  model  development,  surrogate  formulation,  rate  prediction,  etc. 

-  Methods  for  dissemination  and  training 

-  Industry  is  using:  not  predictive  enough  for  computer  design 

•  No  feedback 

•  Little  feedback  even  between  kinetics  and  CFD  communities 

-  Translation  of  uncertainties  in  fundamental  properties/predictions  to 
uncertainties  in  applied  outputs  (engine  performance) 


How  has  the  community  changed? 

•  Today's  perception:  science  will  answer 
everything 

•  1979:  Irvin  Glassman  -  "this  soot  problem  keeps 
coming  around" 

-  Today  -  soot  is  still  a  problem 

•  First  n-heptane  mechanism  ~1970 

-  Still  being  worked  on 

-  First  N-dodecane  mechanism  ~1990s 

•  Similar  fidelity  to  n-heptane 
—  Still  work  need  -  NTC/low-T  regions 

Our  work  is  directed  by  the  problems  facing  industry 


Roadmap/Collaborations 

Common  community  C0  model 

-  Optimized  vs  not 

-  Dynamic  vs  fixed  vs  partially  fixed 

Documentation 

-  Data  must  be  made  available  -  necessary  for  simulation  of 
exp 

-  Models  must  be  available,  documented,  and  consistent 
(e.g.,  forward  and  reverse  rates  consistent  with  thermo) 

-  Should  further  standardization  of  the  requirements  for 
reporting  models  and  experiments  exist? 

Cyber  infrastructure  for  data/models/etc. 

-  Models  and  data  need  to  be  on  web 

•  Improve  the  visibility  of  our  work 

•  Historical  record  of  model  changes 

-  Prime? 


Summary  of  Flame  Chemistry  Discussion  Session 
at  the  1st  Flame  Chemistry  and  Transport  Workshop 


28-29  July  2012 
Warsaw,  Poland 


Chairs: 

Timothy  Ombrello  and  Zheng  Chen 


The  Discussions  Moved  in  the  Direction  of  Trying  to  Understand 
the  Major  Roles  and  Directions  with  Regard  to  Flame  Chemistry 


Validation  of  Flame  Chemistry  Through  Well  Defined  Platforms: 

The  Intermediate  Users  of  the  Chemistry 


What  are  the  Major  Burner  Platforms? 


1.  Flat  Flame/McKenna  Burner 

2.  Counterflow  Burner 

3.  Spherically  Propagating  Flame 


4.  Heat  Flux  Burner 

5.  Micro  Tube/Reactor  Flame 

6.  “Hencken”  Type  Burners 

7.  Others? 


Challenges 

Do  these  platforms  provide  the  detail  and  quantification  needed? 
Can  they  produce  all  of  the  information  that  is  sought? 

-Appropriate  pressures  and  temperatures? 
-1-D,  steady? 

-Optical  and  probe  access? 


Three  Major  Types  of  Platforms  Used  for  Fundamental  Flame 
Studies  —  Their  Strengths  and  Weaknesses 


1.  Flat  Flame  Burner 

-P-0.03-1  atm 

-molecular  beam  mass  spectrometry  measurements 
-thickened  flame  zone  for  better  spatial  measurements 
-how  to  clearly  define  B.C.’s,  especially  when  using  probes 
-losses  to  burner  surface  (heat  and  radicals) 

-higher  pressures? 

2.  Counterflow  Burner 

-P-0.5-5  (10?)  atm 

-ignition,  extinction,  flame  speed  measurements 
-good  optical  access 

-B.C.  problems  under  certain  conditions 
-speciation  measurements? 

-higher  pressures? 

3.  Spherically  Propagating  Flame 

-P~0. 5-100  atm 

-flame  speed  measurements 

-limited  to  flame  speed  and  structure  measurements 


Other  Relevant  Areas  of  Discussion 


Diagnostics 

Quantification  of  new  diagnostic  techniques 

Intrusive  vs.  non-intrusive  measurements,  especially  at  high  pressure 
-spatial  resolution  problems? 

-time  scale  of  measurement  and  quenching  dependence 

Transport  and  Thermo-Chemistry  at  High  Pressure 

Challenge  of  transport  parameters  of  large  intermediate  species 
Can  we  pull  upon  rocket  community  (high  P  and  T)? 

Flames  provide  a  means  of  validation  how  can  diagnostics  help? 

Turbulence? 

What  about  turbulence  and  flame  chemistry  interaction? 
Should  this  be  tackled  at  this  point? 


Defining  Questions  of  Discussion 


1.  Do  the  current  platforms  provide  data  under  relevant 
conditions  (high  pressure  with  and  without  turbulence)  with 
high  enough  fidelity? 

2.  Can  we  think  about  developing  new  platforms  differently  than 
what  was  thought  of  20-50  years  ago  because  of  the 
advancement  of  diagnostic  techniques? 

We  Find  That  There  is  Not  a  Lack  of  Data,  But  Rather  a  Lack 
of  Clearly  Defined  and  Quantitative  Experimental  Data 


How  Can  We  Change? 


The  Way  Forward  -  Two  Paths: 


1.  Extend  current  techniques  for  higher  pressure  with  and  without  turbulence 
and  careful  attention  paid  to  quantification,  systematic  and  measurement 
uncertainties,  application  of  advanced  diagnostics,  etc. 

2.  Develop  new  “standardized”  burner  systems  for  high  pressure  flame  studies 

-similar  to  how  Sandia  developed  Turbulent  Non-premixed  Flames 
(TNF)  library 

-clearly  define  boundary  conditions 

-different  types  of  burners  to  cover  range  of  measurements 
-start  at  1  atm  and  refine  for  higher  pressures 

-heating  and/ or  dilution  at  high  pressure  to  thicken  flame  for  better 
spatial  resolution? 

-some  intermediate  between  a  flow  reactor  and  flame  regime? 

-1-D  and  steady 


The  combustion  and  diagnostics  fields  have  advanced  to  a  point  of 
allowing  us  to  “reinvent”  how  we  think  about  experiments  and  their 
use  for  development  and  validation  of  combustion  chemistry. 


Poster  presentations  at  the  1st  International  Flame  Chemistry  Workshop 


1.  “One-dimensional  (1-D)  Nanostructured  Metal  Oxides  for  Catalytic  Oxidation  of 
Hydrocarbons”  by  Xiaolin  Zheng,  Stanford  University 

2.  "A  database  of  small  molecule  thermochemistry  for  combustion"  by  C.  Franklin 
Goldsmith,  Gregory  R.  Magoon,  William  H.  Green,  Massachusetts  Institute  of 
Technology 

3.  “Further  Developments  in  Multi-Scale  Combustion  Models”  by  M.P.  Burke,  S.J. 
Klippenstein,  L.B.  Harding,  C.F.  Goldsmith,  A.W.  Jasper,  L.  Sheps,  O.  Welz,  J.  Zador, 

H.  Huang,  C.A.  Taatjes,  Argonne  National  Lab.  and  Sandia  National  Lab. 

4.  “Pressure  dependence  in  the  competitive  thermal  isomerization  decomposition  of  the 
cyclohexyl  radical”  by  Enoch  Dames  and  Hai  Wang,  University  of  Southern  California 

5.  “A  Chemical  Kinetic  Study  of  the  Alternative  Transportation  Fuel,  tertiary -BwianoY',  J. 

K.  Lefkowitz  ,  J.S.  Heyne,  Y.  Ju,  F.L.  Dryer ,  Princeton  University 

6.  “Direct  Numerical  Simulations  of  Multistage  Ignition  for  a  Stratified  HCCI  Engine  with 
Exhaust  Gas  Recirculation”,  by  Hossam  A.  El-Asrag*and  Yiguang  Ju,  Princeton 
University, 

7.  “Chemical  Kinetic  Uncertainty  Minimization  through  Laminar  Flame  Speed 
Measurements,”  by  Okjoo  Park,  Peter  S.  Veloo,  Hai  Wang,  Fokion  N.  Egolfopoulos, 
University  of  Southern  California 

8.  A  Comparative  Kinetic  Study  of  Methyl  and  Ethyl  Ester  Oxidation  in  Diffusion  Flames 
by  Pascal  Dievart,  Jing  Gong,  Sang  Hee  Won,  Yiguang  Ju,  Princeton  University 


